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Abstract 

This is the third in a series of works devoted to spectral asymptotics for non-selfadjoint 
perturbations of selfadjoint fe-pseudodifferential operators in dimension 2, having a periodic 
classical flow. Assuming that the strength e of the perturbation is in the range /i 2 <e< h 1 ^ 2 
(and may sometimes reach even smaller values), we get an asymptotic description of the 
eigenvalues in rectangles [-1/C, 1/C] + ie[Fo — 1/C, Fq + 1/C], C^> 1, when eFo is a saddle 
point value of the flow average of the leading perturbation. 
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1 Introduction 

This work is the third in a series devoted to non-self adjoint perturbations of selfadjoint semiclassical 
pseudodifferential operators in two dimensions, whose classical bicharacteristic flow is periodic 
on each energy surface. The previous works in this series are |12l 118) . and more recently, in 
collaboration with S. Vu Ngoc, the authors have begun a study of the case when the classical flow 
of the unperturbed operator is no longer periodic but rather possesses invariant Lagrangian tori 
with a Diophantine property — see |14j for the first work in this direction. 

In this work, we continue with the perturbed periodic case. After switching on a perturbation 
of size e, the spectrum will be confined to a band of width 0(e), and the more precise distribution 
of eigenvalues is very much governed by the flow average of the imaginary part of the leading 
symbol of the perturbation. In the previous works, we studied the eigenvalues associated to non- 
critical values of this flow average or to non-degenerate maxima or minima in a suitable sense 
(after restriction to the 2-dimensional manifold of trajectories in an energy surface). In this paper 
we study the remaining generic case, namely that of eigenvalues associated with a non-degenerate 
saddle point. 

We will work under the general assumptions of ^21, EHj that we now recall. Let M denote R 2 
or a compact real-analytic manifold of dimension 2. We shall let M c stand for a complexification 
of M, so that M c = C 4 in the case when M — R 2 . 

When M = R 2 , let 

P e = P(x,hD x ,e;h) (1.1) 

be the Weyl quantization on R 2 of a symbol P(x, £, e; h) depending smoothly on e € neigh (0, R) 
with values in the space of holomorphic functions of [x, £) in a tubular neighborhood of R 4 in C 4 , 
with 

\P(x,Z,e;h)\<O(l)m(Re(x,0) (1.2) 

there. Here m is assumed to be an order function on R 4 , in the sense that m > and for some 
Co, N > 0, 

m{X) <C (X -Y) N °m(Y) 7 X,Y eR\ (X - Y) := (1 + \X - Y | 2 )i (1.3) 
We also assume that 

m>l, (1.4) 
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and 

oo 

P(z,£,e;/0~ Otf> h^O, (1.5) 

3=0 

in the space of such functions. We make the ellipticity assumption 

\poAx,0\>^m(Re{x,0), \(x,£)\>C, (1.6) 

for some C > 0. 

When M is a compact manifold, we let 

P £ = ]T a a , e (x;h)(hD x ) a , (1.7) 

|a|<m 

be a differential operator on M, such that for every choice of local coordinates, centered at some 
point of M, a Qie (x; h) is a smooth function of e with values in the space of bounded holomorphic 
functions in a complex neighborhood of x = 0. We further assume that 

oo 

a a ^(x; h) <~ a a ^j{x)h? , h — > 0, (1.8) 
in the space of such functions. The semi-classical principal symbol in this case is given by 

PO.eM =$>a,e,o(aO£ a , (1-9) 

and we make the ellipticity assumption 

|po, e (*>OI > ^(0 m , M G ™, |C| > C, (1.10) 

for some large C > 0. (Here we assume that M has been equipped with some Riemannian metric, 
so that |f | and (f) = (1 + |e| 2 ) 1/2 are well-defined.) 

Sometimes, we write p e for p ,e an d simply p for p 0-0 - Assume 

P e =o is formally selfadjoint. (1-H) 

In the case when M is compact, we let the underlying Hilbert space be L 2 (M , n(dx)) for some 
positive real-analytic density /j,(dx) on M. 

Under these assumptions, P f will have discrete spectrum in some fixed neighborhood of G C, 
when h > 0, e > are sufficiently small, and the spectrum in this region will be contained in a 
band |Imz| < 0(e). 

Assume for simplicity that (with p = p t= o) 

p"'(0)nrM is connected. (1.12) 

Let H p = p'^ ■ J| — p' x ■ J| be the Hamilton field of p. In this work, we will always assume that for 
E £ neigh (0,R): 

The H p -&ow is periodic on p^ 1 {E) n T*M with (1.13) 
period T(i?) > depending analytically on E. 
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Let q = \(-§-J e=0 Pt, so that 

p e =p + ieq + 0(e 2 m), (1.14) 
in the case when M = R 2 , and p t — p + ieq + 0(e 2 (£) m ) in the compact case. Let 

I fT{E)/2 

fa) = W^n / 1 ° ^VtH v dt on p-\E) n T*M. (1.15) 

J l ft J J-T(E)/2 

Notice that p, (g) are in involution, = H p (q) =: {p, (q)}. In \T2\, we saw how to reduce ourselves 
to the case when 

p e = p + ie(q)+0(e 2 ), (1.16) 

near p _1 (0) C\T*M . An easy consequence of this, also remarked upon in ^2]> is that the spectrum 
of P e in {z S C; \TLez\ < 6} is confined to ] — S,5[+ie](Req) m i n fi — o(l), (Req) max .o + o(l)[, when 
<5, e, /i — > 0, where (Re g) m i n ,o = min p -i (o)nT* a/ (R e s) an d similarly for (Re (?) max ,o- We shall mainly 
think about the case when (q) is real-valued but will work under the more general assumption that 

Im (q) is an analytic function of p and Re (q), (1-17) 

in a region of T*M, where \p\ < 1/|C(1)|. 

Let A ,F ={/*£ T*M;p(p) = 0, Re (q)(p) = F a }. Assume 

T(0) is the minimal period for the H p -tlow (1-18) 
at every point of Ao^o and Ao^,, is connected. 

The connectedness assumption is for convenience only and can easily be removed. Then So := 
P _1 (0)/exp (RH p ) is a symplectic 2-dimensional manifold near the image A 0j f 01 ^-o,f - We 
consider Re (q) as an analytic function on neigh (Ao._f , So). Assume 

This function has Fq as critical value and the corresponding (1-19) 
critical point is unique, non-degenerate and of signature 0. 

Then Ao._f is an oo-shaped curve, and (q) is an analytic function in a neighborhood of that 
curve (which is the level-curve of Re (q) corresponding to Fq). 

In the following, we may assume that Fq = for simplicity. In Section|2]we shall construct an e- 
dependent canonical transformation k £ which is an e-perturbation of a real canonical transformation 
kq, with 

^,k : neigh ({r = 0}, (TS 1 )^) X neigh (Jf , , C 2 ^) neigh (A ,o, T*M C ), (1.20) 

such that p o kq — cj(t), (q) o kq — (q)(r, x, £) and 

p e o Ke =g( T )+ie(q)(T,x,0+O(e 2 ), (1.21) 

where also the C(e 2 ) is independent of t. Here -Ko,o C R 2 is an oo-shaped curve with the 
self-crossing at (0,0) and (0,0) is the saddle point for the function (a;, £) i— > Re (g)(0, x, £) with 
(q)(0,x,0=F (=0). 

In the present work it seems quite essential to assume that 

The subprincipal symbol of P e= o vanishes. (1-22) 
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(In |12| this assumption was an optional one that permitted to get improved results.) 

After further reductions for the lower order symbols, described in Section 2, we get a microlocal 
reduction of P e near Eo.o to an operator P e (hD t , x, hD x ; h) with symbol 

P € (r,x,e,h) = g(r)+ie^{q)(T,x,0+O(e) + ^ P2 (T,x,0 + jPi + -^ (1-23) 

h 2 

= g(r) +ieQ(T,x,£,e,—;h). 



The operator P e is only microlocally defined near {(t, r, x, £) £ T'S* 1 x T*R; t = 0, (x, £) e Kq } 



but that allows us to define asymptotically its eigenvalues in a rectangle ] — A, -k[+ie]-^, i[ and 



they are of the form 

.9 - ^ - ^) + k e Z (!- 24 ) 

where Wj tk are the eigenvalues near of Q(hk — f 2 - — ^h, x, hD x , e,—;h) in the microlocal space 
Lg, (R) defined with Floquet conditions along the two loops of ifo.o as in J2j- Here 0' = (9±, O2) S 
R 2 with 0j = ^ + fcj G Z, and (So, Si, S2) appear as action differences when quantizing /to, 
while ko,ki, &2 are Maslov indices. 

For r S neigh (0, R), let R(t) be the real analytic curve formed by the values of (g)(r, •). That 
R{t) is a curve follows from (|1.17|) and we see that R(t) is of the form 

Im u> = r(r, Re w), (1-25) 

where r is analytic in a neighborhood of 0. Also, let Pc( r ) € neigh ((0, 0), R 2 ) be the critical point 
of p 1 ► (q)(T,p) (with p"(0) = (0,0)), and let p c (r) € neigh ((0, 0), C 2 ) be the critical point of the 
principal symbol 

^ Q°(r, x, t e, — ) = (-7) + 0(e) + -p 2 (r, x, £) (1.26) 
e «e 

of Q(r,x,S„e,h 2 /e;h) appearing in (JOSJ. Clearly, p c (r) = p°(r) + C(e + Put w c (r) = 

(o)(t, p°(r)) and introduce the exceptional boxes 

B(r) = L; |Re(«;-ti;°(r)| <C (e+-), |Imig - r(r,Rew)\ < £Ml±M l, (1.27) 
i e |ln(e + ^)|J 

for some fixed sufficiently large Co > 0. 

Theorem 1.1 We make the assumptions above and especially Ijl.llfl . i|1.13fl . (|1-17|) . (|1.18|) . i|1.19[l . 

I|1.22|l . and put r k — hk - ^ - ^ , k £ Z. Assume furthermore that h 2 < e < h 1/2 . For C > 
sufficiently large, the eigenvalues of P t in ] — -=5, Aj+ie]Fo — ^,-Fo + ^-[ are 0/ £/ie /orm <|1.24fl . 
where the following can be said about the Wj^: 

The number ofw j>k in B(r k ) is + ^)|ln(e + ^)|. 

Ifuij,k B(r k ), then |Re (wj, k ~ w°(rk))\ > C (e + ^-), with C as in (tL27|) . 
There is a bijection b k between the set of these Wj^ k outside B(r k ) and the union of three sets 
of points away from B(r k ): E cxt = E e> E\ citiat = En, E rightint = E ri such that 

/ -JRe (™— mg(T fc ))| 

bk(w) - w = O I e ch 



In I Re (to -tog (7*) 
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Here E e is a subset of {Re (w — w^Tk)) < — Co(e + — )} and En, E r i are subsets of {Re (w — 

w c( T k)) > C-o( e + ~)} ( or w * ce versa but we only stick to the first option for simplicity) that can 
be described by Bohr-Soraraerfeld conditions 

h 2 1 

b e (w, e, — , r k ; h) = 2ir(j + -)h, j € Z, 9 = e, li, ri, (1.28) 

where 

h 2 °° h 2 

b e {w, e, —, r k ; h)~Y, b e(w, e, — ,r k ; h)h v , (1.29) 

in the space of bounded functions of w,e,h 2 /e,T, that are smooth near (0,0,0) in (e, h 2 /e,r) and 
holomorphic in w for 

1 h 2 
\lm(w-w c (T))\ < -\R e (w-w c (T))\, ±Re (w - w c (r k )) > C (e+ — ), 

with a "— " when 9 = e and "+ " for 9 = li, ri. 
Further, 

6@ is holomorphic in a full neighborhood of w = w c (t), (1.30) 

b» e = 0(\w-w c (T)\ 1 -n, ^>2, (1.31) 

h 2 

b e (w,e, — ,r k ) - 2/iln(-/i), ^-//ln/i, 6 rt -/iln// (1-32) 
are holomorphic in a neighborhood of w = w c (t). 

Here fi is a renormalized spectral parameter defined by w = K t h 2 / e (r k , /i; h), with K given in Propo- 
sitions <|6.2|l , lEU ■ Finally 6q can be described as actions along suitable cycles in the complexified 
cotangent space, see Section ITU1 

Inside the exceptional boxes the eigenvalues Wj. k (for each fixed k) continue to accumulate 
to roughly at most 5 curves where three of the curves are the extensions of the curves carrying 
the E e , En, E r i (defined by replacing 2whj <|1.28fl by a continuous real parameter) and one of 
the new curves, which exists under certain conditions, can be related to barrier top resonances in 
dimension 1. There are at most 2 and at least one point (if we exclude degenerate cases) where 
three of the curves terminate and form a "Y". Away from those points we may have crossings of 
two of the curves (like for instance the ones carrying En and E r i). Away from the "Y" points and 
with some margin, the distribution of eigenvalues can be described by Bohr-Sommerfeld rules as 
in the theorem, and near the "Y" points as well as elsewhere, we can get quite detailed estimates 
for the distribution of eigenvalues. Indeed, the eigenvalues can be identified with zeros of quite 
explicitly given holomorphic functions which in most regions appear as the sum of 4 exponential 
functions, and for such functions it is possible to study the distribution of zeros quite in detail. 
(See Davies 7 for inspiring results in this direction and Hager [U] for quite elaborated results 
obtained in parallel with the present work.) The appearance of Y-shaped eigenvalue distributions 
for non-self adjoint operators in one dimension seems to be quite well known and we refer to 
Shkalikov [23, Redparth ^H] and further references given there, as well as to the recent works by 
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Figure 1: Numerical computation of the eigenvalues of the one-dimensional double well Schrodinger 



4- , ■ ? 
' X + ItX 



operator perturbed by a complex even potential, P t — (hD x ) — : 
h = 0.001 and e = 0.8. When computing the eigenvalues of P e , following |3£ 
operator using the Chebyshev spectral method. 



, in the case when 
we discretized the 



L. Nedelec an d E. Servat and A. Tovbis [20]. The F-shaped eigenvalue distribution is also 
readily observed numerically — see Figure 1 on this page, and also Figure 2 and Figure 3 on pages 
8 and 9, respectively. 

Unfortunately it turned out to be exceedingly difficult to give a concise and precise description 
of what happens inside the exceptional boxes in the form of a theorem in less than several pages, 
so instead we refer the reader to the sections IHHTH1 where this description can be found. 

In Section ED we apply our results to the study of barrier top resonances for potentials of the 
form — x 2 + 0(x 4: ), R 2 3i^0. In the preceding sectional we make a remark that permits to 
improve the domain of validity in the direction of small resonances. This gives an improvement also 
in the applications to barrier tops in |12l ITS*] and allows us in the present work to treat resonances 
E with h 1 - 5 < \E\ < h 1 / 3 for every fixed S > 0, while a direct application of Theorem II . II would 
only give the range ti 1 / 3 <C \E\ <C h 1 ^ 3 . (In this special situation one can say that the lower bound 
e ^> h 2 , can be replaced by e 3> h N ° for every fixed N > 0.) 
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0.08 




Figure 2: Numerical computation of the eigenvalues of the one-dimensional double well Schrodinger 
operator perturbed by a complex odd potential, P e — (hD x ) 2 — x 2 + x + iex 3 , in the case when 
h = 0.001 and e = 0.8. 



When starting this project, we underestimated the amount of ingredients needed, and in order 
to keep the work within a reasonable size, we decided to exclude from the paper the very interesting 
case when there are more than one saddle point on the same connected component of (Re (q) )~ 1 (Fq) 
at real energy 0. The most important case here is probably the one with 2 saddle points arising 
because of an anti-symplectic involution (typically (x, £) <— > (x, — £) in the Schrodinger case). We 
hope to take up at least the 2 saddle point case in a future work (having settled essentially all 
heavy technicalities in the present work). We might then also include the interesting case when 
Hl.lSfl breaks down at isolated points, leading to orbifolds. See Colin de Verdiere-Vu Ngoc 0] in 
the selfadjoint case. 

We also ran into a somewhat unexpected difficulty Indeed, for the one-dimensional operators 
Q, we have a pseudospectral phenomenon leading to an exponential growth of the resolvent norm 
in important regions near the spectrum of these operators. This makes it very important to keep 
the errors in the reduction to the operator P t in 1)1. 23|) exponentially small, so that the accumulated 
error in the global resolvent constructions remains controlled. In |12l 118) we avoided that problem 
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0.08 




x 2 + x 4 



Figure 3: Numerical computation of the eigenvalues of the operator P e = (hD x ) 2 
ie{x 2 +Sx), in the case when h = 0.001, e = 0.8, and S — 0.12. Here S is chosen so that the 
perturbation x 2 + Sx is of the same sign in both of the potential wells, although of different orders 
of magnitude. 



by working in naturally adapted norms where the pseudospectral problems disappeared, but that 
does not seem equally easy to do here. This refined reduction is carried out in Section |21 using 
quasi-norms from the theory of analytic pseudodifferential operators originally due to Boutet de 
Monvel-Kree in the simplified variant of The price to pay is the apparent necessity to 
impose the condition l|1.22f> and the upper bound e -C h 1 ' 2 , that should be compared to the bound 
e < h 5 for every fixed 8 > in OEj or even e < 1 in [23 Ell- 



Acknowledgment. We thank S. Fujiie and T. Ramond for interesting discussions on the rela- 
tionship with the complex WKB-method. Those discussion may be important for future general- 
izations. We also thank the MSRI, Berkeley, for the possibility to advance on this project during 
the Spring of 2003. The first author gratefully acknowledges the partial support provided by the 
National Science Foundation under grant DMS-0304970. 
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2 Reduction to a one-dimensional pseudodifferential oper- 
ator 



Let H C p 1 (0) be a hypersurface which is transversal to the -Hp-directions and such that H can 
be identified with neigh (A 0j o, S ). We can then identify A ,o with a curve H ,o m Ho- 

Let / = g~ x op, where g is the unique increasing analytic function with <?(0) = such that the 
Hf-&ow is 2-7T periodic with 2ir as its minimal period on Ao^- 

Let a : neigh (-Ko,Oj R 2 ) — * neigh (Ho t o, Hq) be a real-analytic canonical transformation, where 
Kofl is an oo-shaped curve, as in the introduction. The existence of such a map with a suitable 
K follows from Theorem 5 in [5] (see also ^2\), according to which if and il 2 are two compact 
real-analytic symplectic manifolds of dimension 2, possibly with boundary, which have the same 
area and such that there exists an orientation preserving analytic diffeomorphism between them, 
then Hi can be mapped onto O2 by an analytic canonical transformation. 

We extend a to a canonical transformation 

K : neigh ({r = 0}, T*S% T ) x neigh (K ,o, R£, C ) -> neigh (A ,o, T*M), 

with / o k = t in the following way: Extend Ho to an analytic hypersurface H in the full phase 
space which intersects p _1 (0) transversally along H$. Let t be the grad-periodic function defined 
near Ao,o which solves 

H f t = 1, with t = on H. 

Because of the 27r-periodicity of the -ff/-flow, we see that t is well-defined up to an integer multiple 
of 2tt. Using that Hf and Hj commute, we notice that if p is a point close to Af^o, then we can 
write 

p = exp (tHf + tH~)(pq), po S H , 

where r € R is small and t £ R is well-defined modulo a multiple of 2tt. 

Put n(t, t; x,£) = exp {tHf + rHj)(a(x, £)). Then n has the desired properties. 

As in Section 2 of 13 j7 we introduce the triple S = (So, Si, S3) £ R 3 of action differences, with 
Sq corresponding to a closed i? p -orbit C p~ 1 (0), and Si, S2 corresponding to the left and right 
closed orbits of the oo-shaped set H ,o- Let 9 = (0q, 61,62), 6j = Sj/(2irh) + kj/4, where kj £ Z 
is a suitable Maslov index. Let L^S 1 x R) be the space of microlocally defined functions u(t,x) 
in neigh ({r = 0}, (T*S 1 )t iT ) x neigh (Ko,o, R 2 ) that are multivalued but 6>-Floquet periodic as in 
|12| . Section 6, or as in J^j. Let U : L^S 1 x R) — > L 2 (M) be a microlocally defined unitary 
Fourier integral operator as in the cited works. 

Repeating the argument in the beginning of Section 3 of |12| . we may and will assume from 
now on that the leading perturbation q in l|1.14|) has already been averaged along the H p -&ow so 
that the leading symbol of P e becomes 

p t =p + ie(q)+0(e 2 ). (2.1) 

The operator 

P £ := XJ- x P t U, 

has the principal symbol 

p e =g(T)+ie(q)(T,x,0 + O(e 2 ). (2.2) 
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At this stage, we get a complete analogue to the situation in Section 3 in ^21- Proposition 3.2 
there extends, and we get a reduction of P e to an operator P e = P e (hDt, x, hD x ; h), also acting on 
Lg, with a complete symbol independent of t, and the principal symbol still of the form 

p e = g(r)+ie(q)(r,x,0+O(e 2 ), (2.3) 

now completely independent of t. 

As in Section 5 of |12) . we use Fourier series expansions in the t- variable and get a reduction 
of (|2.3|) to the family of operators: 

P e (h(k-^)-^,x,hD x ;hj, keZ, (2.4) 

where fco € Z is a fixed Maslov index and it is understood that we only consider such values of k 
for which the first argument (r) in P e is small. 

a) In the general case, without any assumption on the subprincipal symbol of P e =o, we write the 
full symbol of P e as 



P e (T,x,£;h) =g{r)+e 



i{q) (t, x, + 0(e) + -pi{r, x, f) + h-p 2 (T, x, £) + .. 
e e 



(2.5) 



and consider h/e as an additional small parameter. 

b) When the subprincipal symbol of P e =o vanishes, we have the same fact for P e by the improved 
Egorov property of U — see Section 2 of ^2]- Thus (px) e=0 = and we can write 

^Pi(r,x,£,e) = hpi(T,x,£,e). 



Instead of 1(2.5(1 . we get 

h 2 



P £ (T,x,Z;h)=g(T) + e 



i(q}(T, x, f) + 0(e) + — p 2 (r, x, + hpi + h 2 p 2 + 



(2.6) 



depending on the small parameters e, h 2 /e. 
Summing up the discussion so far, we have 

Proposition 2.1 Let P e be as above satisfying the assumptions l|1.2fl . (|1.3|) . (|1.4|l . I|1.5I) . 1)1.60 . 
I|1.10|l . (11.11(1 . (II . 12|) . (|1.13|l . I|1.17|l . (|1.18l) . I|1.19l) . Also assume that < e < h s for some fixed 
S > 0. Then there exist Gq(x 1 ^) holomorphic in some fixed neighborhood o/p _1 (0), an elliptic 
Fourier integral operator U of order 0, with the associated canonical transformation k as above, 
and an h-pseudodifferential operator A — A(t, hD t , x, hD x ; h) of order with principal symbol 
0(e 2 ), such that the operator 

P e = ei A U- 1 e-* G « P <i e^ G oUe-T: A = Ad M , « Gn P £ (2.7) 

has a symbol P e (T,x,t;,e;h) independent of t , modulo O (h°°) . Here Gq = G^{x : hD x ). 

In the general case, we have 1(2.5(1 . provided that h e < h s , and when the subprincipal symbol 
of P t= a vanishes, we have 1(2. ft (I provided that /i 2 <e</i < . 

In this proposition all symbols and phase functions are holomorphic in fixed h, e-independent 
domains. The weight Go in 112.7(1 is used to get a first reduction of the principal symbol to the 
form 1(2.1(1 — see also (18.63(1 in section 8. 
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3 Exponential decoupling 



Since some solution operators to the localized 1-dimensional problems later on will have some 
exponential growth, the decoupling result of Section |2 has to be sharpened in the sense that we 
get some exponential smallness control over the remainders. 

First we need to recall some notions about classical analytic symbols and their associated quasi- 
norms. That was introduced in the pioneering work by L. Boutet de Monvel and P. Kree 0, but 
here we shall use the simplified quasi- norms of |23|. If f2 c is open, a classical analytic symbol 
of order is given by the formal asymptotic expansion, 

oo 

where a k are holomorphic in fl and satisfy the growth condition, 

WK CC 0, 3C = C K > 0; |a*(x,OI < C k+1 k k 7 (z,f) e K. (3.2) 
To such an a, we associate the formal differential operator of infinite order, 

oo 

A(x, C, D x - h) = a(x, £ + hD x ; h) ~ £ h k A k (x, £, D x ), (3.3) 

k=0 

where 

* — ' a\ s 

i+\a\=k 

Let fl t CC f2, to < t < t\ be an increasing family of open subsets with to <ti, such that 

distoo(f] s , d&lt) >t — s, for to < s < t < t\, 
with the distance associated to the £°°-norm. Let fj(Aj) > be the smallest constant such that 

where || • \\ Si t is the operator norm from the space of bounded holomorphic functions on Q t to the 
same space on Q s . Then 

oo 

WWlp-YtftfMj) 



is finite for p > small enough, and conversely, the finiteness of |||a||| p for some fixed p > implies 
that a = ^4(1) is a classical analytic symbol on Q tl - If o-,b are classical analytic symbols on Q 
and we let a{x 1 hD; h) and b(x, hD; h) denote the associated /i-pseudodifferential operators for the 
classical quantization, then the composition of these two operators has the symbol 

oo 1 

a#6~£>* ]T -a?a(x,Z;h)DZb(x,t;h). 

k=0 \a\=k 
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Using that the differential operators A, B compose correspondingly, it was shown very simply in 
|33] that 

|||o#6||| p < |||a||| p |||&||U (3.4) 

implying that the composed symbol is also a classical analytic symbol. 

If we prefer to work with the Weyl quantization, then the same result and proof as in |2.'i| remain 
valid, provided that wc modify the choice of the associated infinite order differential operator to 



(3.5) 



A(x,^,D x>i ;h)u = ^ y(^ a ( D x,^, Dy,^ a(x,(;;h)u(y,r}) 
so that 

A(x, t D Xti ; h) = a(x - -£> 6 £ + -D x ; h). 
In the following, we will consider symbols depending on additional parameters; 

a(x, £, e, h;h) = h k a k (x, £, e, h), (3.6) 

including h (which is then viewed as an independent parameter) and consequently the admissible 
values of p will depend on these parameters. When defining |||a||L in the case when h is among the 
additional parameters, we have in mind some specific representation of the form l|3.6[l . To have 
p tending to as some power of h — and that is what we will encounter, means roughly that we 
consider Gevrey symbols. 

Proposition 3.1 Let l(x,£) be affine and real and let a(x, £; h) be an analytic symbol of order 0. 
Then 

\\\[£(x, hD), a(x, hD; h)]\\\ p < 2p\\Ve\\ oc \\\a\\ p . 

Here and in the following we shall consider [£, a] as an /i-pseudodifferential operator of order 0. 
Proof: The symbol of [£(x, hD), a(x, hD; h)} is equal to 

h h 

-{£, a\ = -v{o x ^)a, 

where v = Hi is the Hamilton field of £. To a we associate the infinite order differential operator 
A(x,^, D x \ h) as in l|3.3|) . Similarly, we have 

OO j oo 

[l(x,hD),a(x,hD;h)]^J2 hjB i( x 'Z> D ^ = jS*** '*' 

3=1 1 k=0 

with Ck — v(d x ^)(Ak) in the sense that v acts as a differential operator on the coefficients of A k . 
Thus, 

Bj = \Ci-x = -Md^iA^). (3.7) 

We can also express this as 
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where t tv denotes translation in the complex domain by the vector rv and r may be complex. 

Assume for simplicity that \v\ = |H|oo < 1- Then for [x, £) € fi s , t < s + \r\ < t < t\ — \r\ and 
for u holomorphic in a larger domain, we get 



u)(x,£)| 



|(Aj_iT r „u)((a;,^)+n/)| < ~ — sup|r ri/ w| 

(t - (s + |r|))a-i nj 

< h-Mi-r) a h ^['L-x Sup 
(t- (s+\r\)y 1 f2j +|r| 

For to < s < t < ti, choose t = t — \r\, 2\r\ < t — s: 

|r_ ri/ o A,_ a oT ru u(x,£)\ < f^A^i)- sup|u|. 

In other words, 

and from the Cauchy inequality and (|3.8|) . we get 



Here we choose 6 so that := minimizes gnj^u-i ■ We find 



BjUt ~ ~6(t-s-2S)>-i [ t _ 8 y _M_ {1 _js_y^ 0<26<t-s. 



Hence, 



so 

/ J -(B i )<2/ J -_ 1 (A J -_ 1 ), 

and 

oo 

PIL<2^/ J -_ 1 (^_ 1 )^ = 2p|||A||| p . 
j'=i 



□ 



In the following, we allow a finite number of families )tj;<t<t£> CC O, f = 1, ...,N as 
above. If |||o|||p denotes the p-quasi-norm, defined with the help of the v-th family, we define 

INIIp^EIIHIIp- ( 3 - 9 ) 

V 

Here it is understood that if a is parameter dependent with h among the parameters, then we use 
the same representation l|3.6|l when defining each of the quasi-norms |||a|||p- Notice that we still 
have l|3.4(l . 
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Proposition 3.2 Let g(x,£; h) be a classical analytic symbol of order 0, defined in a finite union 
D = {J^ =l D u of polydiscs with Q, CC D, and let g(x, hD; h) be the corresponding h-pseudodifferen- 
tial operator. Let ||| • ||| p be a quasi-norm of the form (|3.9[1 with the corresponding family Vl" CC D v . 
Then for p small enough, we have 

\\\\g(x,hD;h),a(x,hD;h)}\\\ p < C(g)p\\\a\\\ p . (3.10) 



Proof: Write g = go(x,0 + hgi(x,t;;h), where gi is a classical analytic symbol of order 0. We 
notice that 

\\\hgi\\\ p <Cp\\\ gi \\\ p , 
(where hg\ is viewed as a symbol of order 0), so on the operator level, we have 

|||[/j<?i,a]||| p < 2|||/i 5l ||| p |H|| p < 2Cp||| 5 i||| p |||a||| p . (3.11) 

Hence it only remains to treat the contribution to the commutator from go(x, hD). We may assume 
we work in a polydisc centered at (0,0) with the radii n, ri, .., r n , si, .., s n . Then, 

go(x,0= E 

a,/3GN" 



where 



'P^s 13 < oo. 



E \9o" 

Now choose the classical quantization for simplicity. On the operator level, 

go(x,hD)=Y,9Z'^ a (hDf, (3.12) 

where the sum converges in the space of analytic symbols, since Hla^lHp < rj, \\\hD Xj \\\ p < Sj+p =: 
Jj , and we can allow some shrinking in rj , Sj and choose p > small enough. 

Using that \\\[xj, a]\\\ p , \\\[hD x .,a]\\\ p < 2p|||a||| p , in view of Proposition 3.1, we see that 

\\\[x a (hD x f,a]\\\ p < {a 1 r a - e ^ + a 2 r a -^ + .. + a n r a -^+l3 1 r a ^~^^..+l3 n ^ 

which can be written more briefly as 

\\\[x a (hD x f,a]\\\ p < (d ri +.. + d rn +d~ Sl +.. + chJ(r a s l3 )2p\\\a\\\ p . 

Hence 

\\\[g (x,hD),a\\\ p < [(d ri +.. + d rn + .. + a 5 „)(El5o' /3 k Q ^)]2p|||a||| p . 

a. 13 

□ 

Remark. The proposition remains valid for the Weyl quantization. Indeed, only l|3.12|l has to be 
modified, by adding a term hgi(x, hD; h) to the right hand side, where g\ is an analytic symbol of 
order 0. 

Returning to the considerations of Section 2, let us consider the analytic symbol, 

P = 9(r) + eq(r, x, £) + hr(t, r, x, £, e; h) = g{r) + eq, (3.13) 
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defined in neigh (t S S 1 , t = 0; ((S 1 + iR) X C) X f2), where f2 C C 2 , ^ is open. Here we assume 
either that r is a classical analytic symbol of order or simply that r has an asymptotic expansion 
in integral powers of h in the space of holomorphic functions. We have already seen in Section 
2 that after a finite number of conjugations, we may assume that r is independent of t modulo 
0{h N ). We also assume that g' ^ 0. 



In the general case, we have 



0(1+*), (3-14) 



and when the subprincipal symbol of P €= o vanishes, we have hr = 0(h 2 + eh + e 2 ), 

q = 0(l + h 2 /e + h + e). (3.15) 
In the two cases, we shall assume respectively that 

h < e < h 5 , (3.16) 

and 

h 2 < e < h 5 , (3.17) 

for some 5 > 0. Then in both cases, we have <f = 0(1). 

We shall see how to eliminate the independence by conjugation with a pseudodifferential op- 
erator up to an exponentially decaying error. The problem can be attacked directly, but it seems 
that we get better remainder estimates if we first reduce ourselves to the case, when 

g(r) = r. (3.18) 

This is possible by means of a holomorphic functional calculus. Indeed, let / = g^ 1 be the 
inverse of the map g. It is easy to see that f(P) is well-defined in the sense of formal analytic h- 
pseudodifferential operators or equivalently in the sense of composition of classical analytic symbols. 
(When r is merely assumed to have an asymptotic expansion in powers of h, we consider those 
hs as additional parameters.) We also see that f(P) as a symbol has the same properties as 
P above, but now with g given by l|3.18|) . It will also be easy to return to the original P, for 
if Ad^P = e A Pe~ A , then at least formally, AdAf(P) = /(Ad^J 3 ), and to say that a symbol is 
independent of t is equivalent to saying that the corresponding operator commutes with translations 
in t and this latter property is stable under taking holomorphic functions of the operator. Until 
further notice g will be given by l|3.18|l . 
Using Proposition 13. 21 we get 

[P,A] = h~{g,A}(t,x,hD t , x ;h)+R(A), (3.19) 
i 

R(A) = e[q,A], 

where 

\\\R(A)\\\ p <C(q)ep\\\A\\\ p: (3.20) 

assuming that the ||| • ||| p -quasi-norm is chosen as in Proposition 13. 21 
Consider the map 

4h Ad A (P) =e A Pe~ A , (3.21) 
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where \\\A\\\ p is supposed to be small. At least formally, 

Ad A (P) = e adA (P). (3.22) 

We expand 

Ad A (P)=J2^A) k (P), (3-23) 
o 

and get the expression for the differential 

5A»a,ds A {P) + jr^ E (ad^ad^ad^ (P). (3.24) 

fc=2 ' v+M=fe-i 

An application of Proposition 3.2 shows that the p quasi-norm of the last term can be estimated 

by 

co 1 oo 1 

C ( P ^Efc! E (2|||A||| p r+l|<U|||, = C(P)p(£ — — (2\\\A\\\ p ) k ^)\\\SA\\\ p 

k=2 ' v+[i=k-l k=2 ^ >' 

= C(P)p(e 2 ^-l)\\\SA\\\ p . 

So, if we assume some fixed upper bound on |||A||| p , we can represent the differential of Ij3.21|) as 

8A i ► ad 5A (P) + K(A, 6A), \\\K(A, SA)\\\ p < C(P)p\\\A\\\ p \\\SA\\\ p , (3.25) 

and combining this with l|3.19|l . (|3.20|l . we get the expression for the differential: 

SA ^ -§{g, 8A} + K(A, SA) = -±g'(T)d t 5A + K{A, SA), (3.26) 
\\\K(A,SA)\\\ p <Cp(\\\A\\\ p + e)\\\SA\\\ p . 

Consider the linear problem 

-g'(T)d t A = B-(B), (3.27) 
i 

where B is a classical analytic symbol of order and 

1 f 2n 

(B) r, x,£;h) = — B(t, t, x, £; h)dt. (3.28) 
27r Jo 

It has the solution 

,-2tx 



(3.29) 



A = C(B) = C(B - (P», C(B) = f (f- \)B{t - s, r, x, Qda. 

g'(T) J 2tt 2 

Clearly (with a convenient choice of the families f2 s ), 

\\\C(B)\\\ p <C Q \\\B-(B)\\\ p <C \\\B\\\ p . (3.30) 

Choose 

P<-. (3.31) 

e 
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We look for A of the form J^c? A? sucn that Ad^(P) is independent of t. We shall do the con- 
struction by successive approximations in a such a way that uniformly during all the steps, 

\\\A\\\ p < 0(1), p\\\A\\\ p < h, and hence \\\K(A, SA)\\\ p < M\\\5A\\\ P , where < 1. (3.32) 

To start with, assume, as we may, that r — (r) = 0(h 2 ) in the sense of ordinary symbols. Choose 

A = C(r), \\\A \\\ P < C \\\r - (r)\\\ p = 0(h 2 ). 

Then, using (|3~2"B|) . |332"|l . 

Ad j4o ( J P) = P-hr + h(r) + hr u 

where P — hr + /i(r) is independent of t and |||?"i|||p < 0\\\r — ( r )lllp- 
Put 

A 1 =C(r 1 ), PJ p <C o 0|||r-(r)||| p . 

Then 

Ad Ao +A 1 (P) = Ad Ao (P) - h{n - (n)) + /w 2 =3 + eg + /i(r) + fc(n) + hr 2 , 

and 

|||r 2 ||| P < e\\\r x \^ <0 2 |||r-(r)||| p . 

Since 9^1 the procedure will converge geometrically and we get a formal solution A with 
Plll^CJIr-^lll,. 

By construction \\\A\\\ p = 0(h 2 ) for < p < min(l, /i/e) and we have defined Ad^-P) by l|3.23|l . 
We define e tA as a formal analytic symbol of order 0, by 

^ t k A k 

fe=0 



so that 

|||e' A ||| P < S A ^% d t e tA = Ae tA = e tA A, e 0A = 1 
in the space of formal symbols. Similarly, we see that 

00 t k xc\ k 

Ft = AdtA = e ^ A= J2t^A 
o M 

satisfies d t F t = ad^ o F t , F = id. We then verify l|3.22|) simply by noticing that 

d t (e tA Pe- tA ) = &d A {e tA Pe- tA ) 
Notice that from the fact that 



p = 0(h 2 ), A = J2 h v A v , a v = A V {1) 



we infer that locally \a v \ < f v {A v ){Cv) v , with £ f u {A v )p» < 0{h 2 ). Hence \a v \ < C h 2 (Cv/p)» , 
so 

\h v a v \ < C a h 2 {CvY(-Y, 

p 
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so A is an analytic symbol with h replaced by h/p and can be realized with an uncertainty 
0{h 2 )e- 1 / ( - ch / p '>'> = 0{h 2 )e- p /( c ^—see [23]. Taking p as large as possible respecting (jS^TJ, 
and recalling that we also assume that p is bounded, we get the uncertainty 

e>(/i 2 )e- 1/c(£+ ' l) . (3.33) 

This discussion can also be applied with A replaced by e A . Let B be such a realization of 
e A . Then we can view A as an analytic symbol of order by declaring that h is an independent 
parameter. Let B~ x be a parametrix, so that Bf^B~ x = 1 + ©(e -1 /^' 1 )), if we also denote by 
B^ 1 a realization. From the construction, it follows that B o P o B^ 1 has a symbol which is 
^independent up to an error of the size l|3.33|l . With this in mind, we can state 

Proposition 3.3 We can construct A in Provosition 12. H such that the symbol P e there is inde- 
pendent oft up to an error which is 0(l)exp (— 1/C(e + h)). Here we assume l|3.16|) in the general 
case and H3.17|l in the case when the subprincipal symbol of P e =o vanishes. 

We end this section with a heuristic discussion explaining why we eventually will assume that 
the subprincipal symbol of P e= o is zero. After decoupling by eliminating the i-dependence as above, 
we get a family of 1-dimensional operators i|2.4[l that we consider at the branching level. If we first 
consider the general case without any assumptions on the subprincipal symbol, we can expect to 
have an estimate on the inverse of these operators or on the associated Grushin problems roughly 
of the order 

-Ks + ^))=- (*(H))- 

In order to combine everything, we would like this quantity times (|3.33|l to be <C 1 . This is clearly 
not the case. 

In the case when the subprincipal symbol of P e= o vanishes we expect to improve the estimate 
on the 1-dimensional resolvents or inverses of Grushin problems to roughly 

xp (c(i + >±!l)) =mq ,(c(i + ±)). (3.35) 

This leads to the condition 

which simplifies to the condition 
which in addition to l|3.17|) gives 



h hi) \ \ h e 



e h 1 
+ - < 



h e e + h 
e 2 < h, h 2 < e, (3.36) 



h 2 < e < h 1/2 . (3.37) 
This is the condition on e stated in Theorem 1.1. 



4 Transition matrix at a branching level. 

In this section and the following one, we study certain model problems. Much of the material is 
standard and close for instance to \TI\, ^B], and [Hj (see also jH] for the C°°-case), but our setup 
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is somewhat different, and we need to recoilect some of the basic facts before returning to our 
operator P e . 
Consider 

{-(xD x +D x x)-a)u = 0, (4.1) 

or equivalently the equation 

d 1 
X-—U = ha )u. 

From Proposition 11 in ^H] we recall that the solutions of 14.1fl in 2?'(R) form a 2-dimensional 

subspace of S'(R). 

For x > 0, we express u as u\x la ~^ , 

For x < 0, we express u as u?,\x\ la ~i , 

For £ > 0, we express u(£) as u 2 £~ M- 2 ; 

For £ < 0, we express u as W4|£|~ m ~2. 

Here u{£) = Fu(£) = ^= J e-"«u(a;)d cc is the Fourier transform and we observe that l|4.1|) is 
equivalent to 

{\{^ + D^) + a)u = Q. (4.2) 
If |Ima| < 1/2, we have two solutions u = U± of (14. 1ft - given by 

= fr(0r* _ia , MO = H(-m^~ ta , (4.3) 

where if = l[o.+oo[ is the Heaviside function, and the general solution to (|4.1|) becomes a linear 
combination, 

u = u 2 U + +UiU-. (4.4) 

We see that [/+ is the boundary value of a holomorphic function in the upper half-plane that we 
also denote by and for x = iy, y > 0, we get 

1 r°° 1 1 

U + (iy) = — / e-y^~ ia d4 = -7=r(i - ia)y ia -i. 
V 27r Jo V 2vr ^ 

Thus for real x, 



which gives 



U+{X) = r( 2" ,a) (— ) ' 

t/+(x) = -i=r(i - <a) x I ef : +<f C7' ? > °. (4.5) 

+ K ' V2^ 2 ; \ e -2 a -' l T\x\ la -2, x < 0. V 7 



Similarly J7_(x) = J7_(x — iO), with 

M"*iO = -±= [ e^|£|---^ = -±=r(± - V > 0, 

V27T J-00 V27T <S 

SO 

U-(x) = -^=r(- - ia)«* - i0)) fa -i, 

V 27T ^ 

M*) = "^ r (o -*«) X 1 _5<*«*l_lia-i ~ ^ n ( 4 ' 6 ) 



/2tt v 2 
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Now let it be a solution of Q4.ip . so that 

u = u l H{x) x -^ +ia + u 3 H(-x)\x\-? +ia = u 2 U+ + U4U-. (4.7) 
Using 14.5|l and (|4.6|l . we get 

ui = ^Lr(i - ia)e% a+i *U2 + -^=T(- - iaje^^^u^ (4.8) 
2 2 



W3 = — !=r(^- - ia)e 5 Q 'J U2 h F=r(i - ia)e*° +1 *M4. 

2 2 



Here we want to express 7/2, it 1 in terms of 113,1*4. From i|4.8[l . we get 

"2 \ = / r(i- lQ ) e 

\ r(5+ia) 
where we also used the reflection identity, 




(4.9) 



Recall that T(z) is meromorphic with simple poles at — k, for k E N, and no other poles. The 
reflection identity above can also be written 

r(z)r(i-z) = -^, 

smirz 

and implies that if T(z) = 0, then 1 — z has to be pole, so 1 — z = —k for some fcgN, z = k + 1, 
which is impossible since we also know that T(k + 1) = k\ ^ 0. Hence T(z) has no zeros, and 
is entire. The transition matrix in (|4.9|) has determinant 1 and is an entire holomorphic function 
of a. The relation (|4.9|l remains valid therefore for all a £ C. 

We next compute the transition matrix analogous to the one in (|4.9() in the semi-classical case, 
for solutions of 

{-{xhD x + hD x x) - fi)u = 0. (4.11) 

This is of course the same equation as (|4.1|) with a = [i/h. We now require 

u = uix l %~^ , x > 0, 
u = u 3 \x\ l %~^ , x < 0, 

wo = U2 r^~K e >o 

where 

T h u{0 = -L= J e-^' h u{x)dx = ^«(|). (4-12) 
A simple computation gives 



ui/ I e <+*f lief / V M4 / V ai ' 3 a W W 

We summarize the discussion above in the following proposition. 



(4.13) 
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Proposition 4.1 Let p 6 C be such that \Imp\ < /i/2. If u E T>'(R) is a solution of l|4.11|l then 
u € 5'(R) and £/iere exzsi «i, M2, "3, «4 € C such that 

u = u 1 H{x)x i t-i + u 2 T^ (ff(C)r^"') = u 3 tf (-x) M** - * + u^ 1 (h(-£) l^ - ** - *) ■ 

Here Th is the semiclassical Fourier transform defined in (|4.12l) and the coefficients u 2 , Ui can 
be expressed in terms of u 3 , ui by (|4.13|) . The transition matrix which occurs in (|4.13|) is entire 
holomorphic in p and has determinant 1. 

We finish this section by the following observation, which will be useful in section 6. The 
operator Pq = h{xhD x + hD x x) has the principal symbol po(x,£) — x£. For \i 6 C, <C 1, 
define pj £ Pq 1 ^) by pi = (l,fx), pi = (p, 1), pz = (-1,—fi), pi = (-fx, -1). Working in the 
semiclassical limit, define the microlocal null solutions ej of Pq — p, for j = 1, ..,4 by e\ = x l h~2 
near pi, e\ = near p 3 , = near pi w Kp h p2, e 2 = near ^4, e 3 (x) = ej_2(-x), 

j = 3,4. Here is the map (x, £) 1— > (£,—2;) associated to Th- Then a general null solution of 
Pq — p can be written either as 1*262 + u\ei or as 1*363 + ^464, where (14.131) holds. 



5 Asymptotics of the transition matrix 

In this section, we shall derive asymptotic formulas for the entries of the transition matrix l|4.13|l . 
We shall use the following version of Stirling's formula (HZ|), 



T(z) 



2tt 



e- z z z -i(l + — 



1 



1 



12z 288z 2 



00 , I arg z I < 7r — 5, 



(5.1) 



for every fixed d > 0. We apply this in two cases: 

Case A. We have \p\/h ^> 1 and p £ a conic neighborhood of the negative imaginary axis. Then 
ijl avoid a conic neighborhood of R_ and we can apply (|5.1|l . to get: 



V2tt V 2 


^) 
/V 






V 2 










exp 


1 

"2 H 


■ M 








exp 


1 

"2 H 


■ M 








exp 


1 

:^ 




case we have 












V2^ r( 2 


— i 




exp 


ip 

~h ~ 


^ ln( 



f rM^-ir) + o(-) 

h A h p 
re /i 2i/x /i 



:-) 



(5.2) 



Here and in what follows In always stands for the principal branch of the logarithm. 
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Case B. We have \fJ-\/h S> 1 and fi avoids some conic neighborhood of the positive imaginary axis. 
Then we can apply the earlier results with /x replaced by —fj, and get: 



=r(i 



exp 



■— + «- ln(t/z) - — In ft, + ©+(-, 
ha a fi 



(5.3) 



If Re > i|Im^|, we combine this with the reflection identity: 



2tt 4 V l 2 + V 2cosh^ 



and the fact that 



ln(i/i) — ln(— ifi) = 



in this region, to conclude that 



(5.4) 



It is straightforward to establish a corresponding estimate in the region Re fj, < — ^\Im fi\, and we 
can summarize both cases in: 



0+{-) + = 0(e- 2 ^ Re ^ h ), |Re/i| > A|Im 



/' 



C 



(5.5) 



Remark. When ji is real we have 



ReO+(-) = ReO_(-) = 0( e - 2,r W/' 1 ). 



(5.6) 



In fact, if we first assume that /i h, we get from (|5.2II . (|5.3|l : 



r(|-^) 



exp 



*r -«7-ln/i + i r ln/i- — + £)_(-) 
h h n In ji 



2 1 "/i' 



exp 



^ + ^ in/i -^ in ^2ft +o+ y. 



and using that the second quantity is equal to the complex conjugate of the other, we conclude 
that (|5.6[1 holds in this case. In the case fj, <C — h, we can use the same argument. In this case 
ftTSjl . JOl give 



r(j-<£) 



exp 



= exp - In | M | - In h + ^ + 



/' 



and we can again conclude that Q5.6|l holds. 
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In any closed sector away from «R_ , we can use (I4.13f) and (|5.2I) to get 

a 2 , 3 = exp f 1 ^ ln(-i M ) - - p + ^ + ^-(^)J J • (5-7) 
In any closed sector away from zR+, we can use (I4.13f) . the reflection identity and l|5.3|) . to get 

a23 = 2cosH^)e^ l ^-^-^-* ho ^ (5.8) 

Using l|5.3|) . I|4.I3(I . we get for /i away from a sector around zR + : 

ai 4 = eH-nHiri-izp+n-'F+ihO+li)) (5.9) 

In a sector Im/i > — C|Rep|, we use the reflection identity 

- 2 71 -2 cosh 7T^-, 



r(§ + i£) h' 

to get 

ai 4 = 2 C o S h(^)e^ u ^- i ^ + ^- ih0 -^. (5.10) 

Combining the asymptotic formulae H5.7|) - (|5.1()J) . we may summarize the discussion in this 
section in the following proposition. 

Proposition 5.1 We have the following tableau for the coefficients of the transition matrix 14.13fl . 
when \fi\/h 3> 1. In all cases: 



For Rep > C 1 |Im/x|: 



For Im /i > — C|Re/z|: 



ForRefi<-C 1 \lmfi\: 



For Im (i < C|Re fi\: 



a2,4 = -e^^ 2 , oi,3 = e* h+l * 



a 2 3 = e i(MlnM-iTM-M+^+*C , -($)) ) 



02i g = e *(fInif-i^-/H-^-HhO-(*)) ) 

a X4 = e i(-Min^- l ^±^ I+A1 -^- 2 ho_(^))^ 



a 2i3 = e *^ to (-/*)-iH-^-HfcO-($)) j 
0l 4 = e U-»M-f*)+»-"r L +ihO+(%)) , 



02,3 = e 



fc(p\n(i l i)-iS£±iwii-li+*£—ihO + (%)) 



ai 4 = e t(-Min(v)-^+M-^+^o+(ir))_ 

Here the terms with iiir/j, in the exponents indicate that we should take the sum of the two 
possible terms with the same remainders 0+ or O- in the exponent for each of the two terms. 
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6 The one-dimensional spectral problem 



We now return to Proposition ^. II which shows (when combined with the exponential remainder 
estimates of Section [3} , that the study of P € near Ao,o (considered in a suitable weighted space) 
can be reduced by conjugation to that of P e acting on the space L^S 1 x R) of functions defined 
microlocally in some fixed neighborhood of r = 0, (x, £) £ -Ko,o C R 2 in r*^ 1 x T*R, with a 9 = 
(#0j$1j$2) Floquet periodicity condition. Here -Ko,o is an oo-shaped curve as in the introduction, 
and 9 was defined in the beginning of Section [21 In order to fix the ideas, we assume that we are 
in the case when the subprincipal symbol of P e= o vanishes, so that the symbol of P e is given by 
l|2.6ll . and we are then in the parameter range h 2 -C e < h s . At least formally, the study of P e can 
be reduced to a family of one-dimensional problems by a Fourier series expansion in the t- variable 
and (as noted in (|2.4[l ) we get the one-dimensional operators 

P e (h(k-9 ),x,hD x ;h), k G Z, (6.1) 

where 9q = 5o/(27r/i) + fco/4, and we restrict the range of k by requiring that h(k — 9q) be small. 
The operators i|6.1|) should be considered as acting on the microlocal space L 2 ,, (R) defined similarly 
to Ll(S l x R), with 9' = (9 1 ,9 2 ). Using lj2~S|). we see that JOJ becomes 

h 2 

g(r) + ieQ{r, x, hD x , e, — ; h), T = h(k 9 ), (6.2) 
e 

where 

h 2 h 2 
Q{t, x, £,e, — ; h) ~ Q (r, a;, £, e, — ) + /iQi + /i 2 Q 2 + ■•• 
e e 

is holomorphic with respect to (t, x, ^) in a fixed complex neighborhood of {0} x iio.o and depends 
smoothly on the other parameters. We further notice that 

Qo(T,»,C,0,0) = (?)(r,i > (6-3) 

is equal to the trajectory average of q, expressed in suitable real symplectic coordinates, and we 
know by construction that 

(q}(0,x,0=0onK , . (6.4) 
We also recall the assumptions (|1.17|) . I|1.19() . which say that 

(q)(T,x,t)=f{T,Re(q){T,x,0), Re/(r,r) = r, /(0,0) = (6.5) 

^ e (<l)(x,$),(x,£)(Q> 0, 0) is non-degenerate of signature 0. (6.6) 

Here we assume for simplicity that (0, 0) € -Ko,o is the branching point. 

In the following we shall discuss the spectrum of the 1-dimensional operator 

h 2 

Q = Q T = Q(T,x,hD x ,e, — ;h). (6.7) 

e 

Since this operator is only defined microlocally and up to an error O c('+ h ) J , we have to keep 

in mind that for the moment the eigenvalues will be defined only formally and up to errors of at 
least the same size. 
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A first localization of the spectrum. Assume first that (q) is real-valued. Then from the sharp 
Garding inequality, we see that the spectrum of the operator 1)6. 7fl in the band |Rez| < 1/0(1) is 
contained in 

ze C; \Rez\ < 1/0(1), \lmz\ < 0(l)(h + e + ^-)J . (6.8) 

Under the more general assumption (|6.5|l . we see that (|6.8|l can be applied to g(r,Q T ), where 
g(r, ■) = / _1 (t, •). So in the general case, we see that the spectrum of the operator l|6.7|l in the 
band |Rez| < 1/0(1) is contained in 



E r := iz G C; \Rez\ < 1/0(1), \z-f(r,Rez)\ < 0{l)(h + e + — ) j . (6.9) 

We also have 

ii (Or - zr'w < di^^) » for i Rez i < vo(i), ^ ^ e t . (6.io) 

Normal forms near the branching points. Let Qq = Qq(t, x, £, e, — ) be the principal symbol of Q. 
Following J21> |10| - we get the following adaptation of Proposition 5.3 of |12| : 

Proposition 6.1 We can find a canonical transformation: (x,£) i— > £) depending ana- 

lytically on t and smoothly on e, /i 2 /e iwt/i values in the holomorphic canonical transformations: 
neigh ((0, 0), C 2 ) — * neigh ((0, 0), C 2 ), and an analytic function k t \ l i / c (t, q) depending smoothly on 
e, h 2 /e, such that 

Kr,c,h°/e{0, 0) - (x(r, e, /i 2 A), e(r, e, h 2 /<0) 
is the unique critical point close to (0,0) o/ (x, £) i— > Qo{t, x, £, e, h 2 /e) and with 

h 2 

Qo(T,K r h 2 /e (x, £),e,— ) = k h 2 (r,xf). (6.11) 
Moreover, k t ,o,o is real when r is real and 

^-Refc h2 ( T ,0) >0. 

After a conjugation by an elliptic Fourier integral operator associated to the canonical trans- 
formation K rj£ .^2/ e we may assume that the leading symbol of Q T , 

Qo = Qo(r,x,C,e,/i 2 A) (6-12) 

is a function of r, e, h 2 /e and x£. We can get a complete normal form by making further conjuga- 
tions by analytic pseudodifferential operators of order in such a way that the complete symbol 
also becomes a function of r, e, h 2 /e and x£. This is carried out in Appendix A. We get the 
following result which is very close to one of the main results of the appendix b in 

Proposition 6.2 We can quantize K; T ,e,/i 2 /e by on elliptic Fourier integral operator U = U TiCi h 2 /e 
with an analytic symbol, depending holomorphically on t and smoothly on e,h 2 /e, such that 

U- 1 QU = K e>h 2 /e (T,I;h) + 0(e-^), I = P = hx o hD x + hD x o x), (6.13) 

where K e ^%i e {r, L\ h) is a classical analytic symbol of order depending holomorphically on t and 
smoothly on e,h 2 /e. The leading part of K is k e h?/e (jt L ) appearing in Provosition Wi.W 
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The quantization condition. We start with a side remark about normalization. When P(x,hD x ) 
is a selfadjoint one-dimensional /i-pseudodifferential operator of principal type, and z £ R, then 
we can normalize microlocally defined solutions of (P — z)u = 0, z € R, by imposing that 

C-[P, X ]u\u) = \. 

here \ = x{ x i is defined near a piece of the real characteristics and has the property that V% is 
of compact support near the characteristics of P — z and x increases from to 1 when we progress 
in the Hamilton flow direction. (See 52 •) ^ is easy to check that if we view u as a solution of 
(f(P) — f(z))u = 0, then we get the corresponding normalization 

(■^[f(P),x\u\u) = g{z,z){-[P,x]u\u) = g{z,z), 

where 

f(P)-f(z) = (P-z)g(P,z). 

If we drop the requirement that P be selfadjoint or just let z become complex, there is no obvious 
normalization of null-solutions of P — z, but we still have a well-defined scsqui-linear form on 
N{P -z)x Af(P* - z), given by 

Cj\P,x]u\v). 

If we have some additional information allowing us to identify the two null-spaces, then this can 
still be used to normalize null-solutions of P — z. In the following we abandon the attempt to 
normalize completely the null-solutions, since already the operator <2 e= o,h 2 /e=o is n °t necessarily 
selfadjoint. 

By Proposition ^ . 2l we have an analytic symbol f(-;h) depending analytically on r and smoothly 
on e, h 2 /e, such that 

U- 1 f(Q;h)U = P Q . (6.14) 

Notice that if u is a null-solution of Po — /i in a full neighborhood of (0, 0), then (Q — z)Uu = 
near (0, 0), where the spectral parameters are related by 

f(z;h) = n. (6.15) 

Recall from the end of Section 0] that Po — /i has the four characteristic points pj, j = 1,2,3,4 and 
that this operator has the microlocal null solutions ej described after Proposition ^. II When /i is 
real, we check that ej is normalized near pj. If v is a global null-solution of Po — p, with v — VjCj 
near pj, then by Proposition 4.1 we have: 



v% \ _ I a 2 , 3 a 2j4 \ / v 3 
Vl J V Oi )3 a 1A / I « 4 



(6.16) 



Assume for simplicity that k = n T ,t,h 2 /t is defined in a suitable domain, containing pj, j — 
1, ..,4. Let a.j — n(pj), fj — Uej. Then if u is null-solution of Q — z near the branching point 
k((0, 0)), equal to Vjfj near a.j, then (|6.16(l still holds. We may wish to renormalize the fj, by 
putting 

fj - e idi 9j- (6-17) 



27 



Then a straightforward calculation shows that if u is a null-solution of Q — z near the branching 
point, and u = Ujgj near aj, then 

uA = fc 2 , 3 c 2A \ fu 3 \ = e -m-^ a . k . ( 6 .i 8) 

Ul ) V C 1.3 c l,4 / \Ui J 

Here is a natural example of such a renormalization. Assume for simplicity that near aj the set 
Qq (z) takes the form £ = \j(x), where \j is analytic and depends analytically on the parameters 
e, h 2 /e, t, z. Then choose gj so that microlocally near aj we have the standard WKB-form: 

9 j = bj{x;h)ei^ x \ (6.19) 

where bj is a classical elliptic analytic symbol of order 0. The function ipj solves the eikonal 
equation 



dip 

— Xj(x) = 0, with the extra condition ipj(n x (aj)) = 0, 



and bj,ipj depend analytically on the additional parameters r, z and smoothly on e.h /e. Using 
an explicit representation of U we write near aj for j = 1,3: 

f 3 ( x ) = h-^- [ [ ei^ {x ^' e)+ ^ {v)) A{x iy ,e)a 3 {y)dyd9 = e&bjfah). (6.20) 



Here the last equality follows from stationary phase, ip is a non-degenerate phase function gener- 
ating k and near pj we write ej in the WKB-form 

e j (y) = a j (y)ei*iW = |#"' = \y\-ie*» u M. 

The function ipj(x) in appears as the critical value in the stationary phase expansion of 

H6.20fl and solves the same eikonal equation as ipj. 
For j — 2, 4, we get near aj: 

fj(x) = [[[ e^ {x - v ' e)+m+,t '^ ]] A{x 1 y,e)a ] {ri)dyd7 ] de = e&bjfah), (6.21) 



where aj and (pj appear when writing Th&j on WKB-form near Kj^(pj). 

In this case we see that dj = dj(h) is a classical analytic symbol of order 0, depending analyt- 
ically on the additional parameters and with the imaginary part of the leading symbol vanishing 
when Imp — 0, e = 0, h 2 /e — 0. The leading part of dj can be further described in terms of 
symplectic geometry. 

Put 

6j, k ■= dj - d fc . (6.22) 

We have the obvious relation 

#2,3 + #1,4 = #1,3 + #2,4- (6.23) 

Now we work in a full neighborhood of Kq^q (introduced after H1.21|) l. Recall that we have the 
points a±, a 2 , on the four crossing branches of Kq q distributed with positive orientation 

around the branching point. We may assume that a 3 , a a are situated close to the left closed curve 
7 1 of Kofi and that ai, a 2 are situated close to the right closed curve 7 2 of ifo.o- Start with a 
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microlocal null-solution to Q — z near k((0, 0)), of the form 1x4^4 near 04 and of the form u 3 g 3 near 
a 3 . Here for the moment 1*3, 1*4 can be prescribed arbitrarily, and we then know that u = Ujgj near 
aj for j = 2, 1, where 112, u\ are given by (|6.18|l . We require temporarily that u is a well-defined 
single- valued null-solution along the whole left closed component 7 1 of -Ko,o- Then if we follow u 
around the exterior part of 7 1 from 124 to 03, we get 

u 3 = ei S3A u 4 , where S 3A = / £dx + O(h) = S$ A + O(h), (6.24) 

with 73 4 denoting the exterior part of 7 1 which joins a 4 to 0:3. 

Now recall that we really want u G Lg,, 9' = (0i, 6*2), where 9j = tJ-^ + ^- where Sj is a real 
action difference related to the reduction in Proposition 12.11 and fcj£Za corresponding Maslov 
index. This means that u should be multivalued, but Floquet periodic along 7 1 in the sense that 

7> = e~ 2 ™V (6.25) 

where j^u denotes the extension of u along one loop of 7 1 which we assume to be oriented in the 
following way: a.4 — > a 3 — > (0, 0) — > a.4. Starting near 03, we get j}u near the same point in two 
steps: 

U393 -> "434 -* e^ 53 ' 4 M433- 

The Floquet condition (|6.25(l therefore becomes e~ 2mSl u 3 = e lS ^ A ^ h u 4 , or equivalently 

u 3 = e 2m ^ + ^ s ^Ui, (6.26) 

instead of (|6.24|) . 

Similarly, let 7 2 be the right hand loop in -Ko.o with the orientation: 0.2 — > ct\ — > (0,0) — > «2- 
Then, if we want u to extend to a null-solution in Lg, near 7 2 , we get the analogue of Ij6.26|l : 



Ui 



= e 2m92+ -k Sl - 2 u 2 , (6.27) 



with Si t 2 defined as in (|6.24() with 73^ there replaced by 71,2, the exterior segment in 7 2 that joins 
a 2 to ct\. 

Start near 014 with 114(74, use (I6.26|) to get u 3 and then (|6.18|l to get 1*2, u±: 

U 2 = (C 2 ,3e 2 ™ 9l + £ S3 < 4 + C2, 4 )"4, , fi 2 o^ 

ui = (ci, 3 e 2 ™ ei + ^+ CM K ' 1 ' ' 

and in order to get a global solution in Lg,, we also need to apply (|6.27() . which gives our global 
one-dimensional quantization condition 

= C2, 3 e 2 ^ +e ^+^ s ^ +s ^ + c 2 , 4 e 2 ™ f,2+ ^ Sl < 2 - c h3 e 2me ^ Ss '* - c M , (6.29) 

where we took U4 = 1. 

In this relation, we substitute (|6.18() . I|6.22[l and get after multiplication with e^ 1 - 4 ''*: 

= a 2 , 3 e^ (§3 ' 4+Sl - 2) + a 2A e^ 1 - 2 - ai, 3 e^ §3 ' 4 - ai, 4 , 
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with 



S h2 = Si, 2 + 6>i,2 + 2irh0 2 =Si, 2 + 0i,2 + S* 2 + hk 2 ~ (6.30) 

•33,4 = ^3 4 + 03 4 + 27l7l0i = $3,4 + 03,4 + Si + hk\ — , 



where we recall that 



With 



Sj_ kj_ 
271- 4 



~ ^ 7T 

Si,* = Sj,Jfc + /l— , (6.31) 



we get 



= a 2:3 e^ (S3 ' 4+Sl ' 2Mf + a 2 , 4 e^ Sl < 2 - a^e^ 3 - 4 - a M e^ . (6.32) 



Proposition 6.3 Assume that 

h 2 h 2 h 2 

Q(r,x,£,e, —;h) ~ Q (r,a;,C,e, — ) + hQi(r,x,£,e, —) + ... 

is holomorphic in (r,(x,£)) € neigh(0, C) x neigh(ifo,Oi C 2 ) and depends smoothly on e, — € 
neigh(0, R). /fere -Ko.o is aw oo-shaped curve with the self-crossing at (0,0). Assume furthermore 
that 

Qo(t, x, i, 0, 0) = (q) (t, x, - /(r, Re (q) (r, x, 0), 

where f is an analytic function with /(0, 0) = 0. We assume next that along K , Re (g)(0, x, £) 
I) and i/iai Re (q)(0, x, £) has a unique critical point on K a o , (0,0), which is a non-degenerate 
saddle point. When z £ neigh(0, C), put 

M = f(z;h) 

where f(z;h) is an analytic symbol introduced in Propositions 6.1 and 6.2, and in ()6.15|) . Then 
z is a quasi- eigenvalue of the operator Q (t, x, hD x , e, hj acting on Lg,(R) if and only if the 



corresponding /i satisfies (|6.32l) . In (I6.32[) . f/ie coefficients ai.3, ai,4, 02, 3, and 02,4 are introduced 
in Proposition 4.1, and i/ie quantities Si, 2 and £3,4 are defined in (|6.30(1 and H6.31|l . TTiey depend 
holomorphically on /i with 

d„Sj, k = 0(1), (6.33) 

and w/ien /i is rea/, we Ziawe 

lmS itk = o(e+!jJ . (6.34) 

In the formulation of the proposition, we leave the notion of a quasi- eigenvalue undefined and 
refer the reader to Section 11 for a complete justification of this terminology. 
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7 Zeros of sums of exponential functions 

Here we elaborate on arguments in [7], and a related and even more general discussion can be 
found in Hager 9 . The results established in this section will be used in Section 10. 

Let 71,72, •••,7jv be compact C 1 segments in C such that jj starts at s 3 _i S D(zj-i, rj-i/2) 
and ends at ej S D(zj,rj/2), where we use the cyclic convention and view the index j as an element 
of Z/NZ. We assume that N € {1, 2, ...} is fixed, but allow jj, Zj,rj, Sj, ej to vary with the semi- 
classical parameter h while all estimates below will be uniform in h. Let / be a holomorphic 
function defined in L)j r S^ 1 (D(Zj, rj) U neigh (73+1)), such that 

f = eiSj{z)+ o(i) on7 . 5 

l/l < e i(-imS,W)+o(i) onD^.rj), 
where S*j is holomorphic in neigh (7^) U D{zj-\ 1 rj-i) U D(zj,rj) and 

Im (Sj+i - 5,-) = on D{z j ,r j ). 

In D(zj,rj) we can write 

/(*) = e* s 'W ffJ (z), |. 9j (z)|<0(1). 

We further know that |5j(ej)| > 1/0(1). Standard arguments (see for instance [Sj), including 
Jensen's formula, imply that the number of zeros of gj in D(zj,rj/2) is 0(1) and if oij is a 
segment in D(zj,rj/2) from ej to Sj which avoids the zeros Wi, ■■■ 1 wm of gj in D(zj,Tj /2) such 
that |vararg Q .(2; — u)fe)| < 27r for every fc, then 



(7.1) 



Re^r/ ^dz = 0(l), 
and consequently 

Rg 2^ J a J dZ = 0il) + 1 khj a ReS '^ dz = °M + ^h Re Pite) " (7-2) 
Let 7 be the closed contour given by 71 U ai U 72 U ... U 7jv U ao- We want to study 

1 f f'(z) 1 
N(fn) ■■= 77- / ~7TT dz = ^vararg 7 (/). 
2tti 7 7 f(z) 2tt 7 

When 7 is the oriented boundary of a bounded domain T, where / is holomorphic, then N(f, 7) is 
the number of zeros of / inside T. 

Along jj, we write / = e* § ^ 2 \ ^-(z) = ^(z) + 0{h). Then, 

2^ /, J dZ = 2^ /, = ^h^ ] ~ 

so 

E^- / 7^ = ^Eto( e ^-^iteO) = ^Etote)-^+i(^))+o(i), 
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and hence in view of (|7.2|) and the uniform boundedness of N: 



N-l 

N(fn) = g-j- £ ( 7 - 3 ) 

Here we recall that Im(<Sj - S^+i) = O(ft-) in D{z J1 r j ). It follows that V(Sj - Sj+i) = 0(h/rj) 
in D(zj,rj/2) and consequently that 

Sj(sj) - S , j+1 (s i+ i) = S^z) - Sj+i(z) + 

for any other point z € D(Zj,rj/2). Thus finally, 

AT-l 

^(/,7) = rrE^W ~ Sj+iK))+0(l), (7.4) 

with iDj G D(zj,Tj/2) chosen arbitrarily. Here we can further replace Sj(wj) — S^+i^iUj-) by its 
real part, since ImSj — Sj+i = in D{Zj,Tj). 

8 Skeleton in the region ^> h. 

We now return to the situation in Section [(J We are interested in the solutions /i of l|6.32|l . In 
the following, we will write Sj t k instead of £j,fc, so we are interested in the zeros of the function 
i*b(/x; h) appearing in 1)6. 32|) . given by 

F (n; h) = e* (Sl ' 2+S3 - 4) "^a 2 , 3 + e* 5l - 2 a 2:4 - ai, 3 e^ 53 < 4 - o M e^. 

Pulling out a factor e~ J7r / 2 , we get the new equivalent function 

F{n\ h) = e^ (Sl - 2+S3 ' 4) a 2 , 3 + e^ 5l - 2+7r £ + e^ S3 - 4+7r £ + a 1A , (8.1) 

which has the same zeros as Fq. Here we have also used the explicit formulae for 112,4, ai.3 in 

(S3) . 

Using the results of Section 5, we shall now look at the asymptotics of -F(/x; h), when ^S> 1. 



Case i: Assume that 



C/i < < 1, 



7T 

arg/i 



1 



2 

(The case 2, given by |arg fi+ -| \ < ir— will be reduced to the case 1 by a symmetry argument.) 
In this region, we have l|5.10[l : 

01,4 = 6°-^ ( e ^(-^ ln f+^T+^) +e i("Mlnf+M-^-i 2 F)^ 5 

and we get using also l|8.1|l and 1)5. 7[) . 

+e O-(^)+*(-Mln*+M-^+i^) +e O-(^) + t(-pln^+M-^-^) 

= e^G{n;h), 
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where 



G(n;h) = a 1 +a 2 + a 3 + a 4 , a 4 = a 4 + + a 4 - (8.3) 
ai = e i(-Si,2+S 3 ,4+^lnf-^+^)-0_(i) 

a 3 = e i s 3A+TK 



a 4 ± 



We have 
where 



= e r '/ h , j = 1,2,3, 4±, (8.4) 



ri : = — Im 5i 2 — Im- S*3 4 + (Im /i) In - — - — Re /i arg — + Ira u — hRe O- ( — ), 

|/x| i n 

7T 

r 2 : = -ImSi i2 + — Re/z, 

7T 

r 3 : = -Im5 3:4 + -Re/x, 

: = — (Im /i) In - — - + Re /x arg — — Ira /i ± 7rRe fi + hRe O- ( — ) . 

I A*l * V 

Notice that a 4 ± is dominating over a 4 T when ±Re /x > 0, and in each half-plane ±Re /i > 0, we may 
associate a 4 to the dominating term, modulo an error which is 0(e- 27r ^ etl ^ h ) times the leading 
term. Also notice that the last equations take the form 

n = (Im fi) In -i- - Im Si 2 ~ Im 5*3,4 - Y(/j), (8.5) 

ImI 

7T 

r 2 = -Ito5^ 2 + -Re/x, 

7T 

r 3 = -Im 63,4 + -Re/x, 

r 4 ± = — (Im/z)ln- — - ± 7rRe /z + Y(li), 
H 



where 



Y{p) = (Re/x)arg(-) -lmfi + hReO-(-). (8.6) 
i /x 



Following the general principles, as explained for example in [7] (see also PQ), we shall now 
look for the curves IV*., j, k = 1,2,3,4 , where \aj\ = |cifc|, and we shall especially be interested in 
those parts of Tjk where |<Zj| = \ak\ is dominating over the other \a v \. In doing so, let us remark 
first that we will not see any zeros of G generated by the zeros of a 2 + 03 as a dominating part of 
G, for when Re xz > 0, then r 4 dominates over r 4 and r\ + r\ = r 2 + r 3 and clearly we cannot 
have r 2 = r 3 > max(ri,r 4 f ) + Const. Now when Re/x < 0, r 4 dominates over r 4 and 

?*2 + r 3 - 7rRe /x = r x + r 4 , 

so that r 2 + r 3 < r± + in this case, leading to an even stronger conclusion. 
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We now begin look at the location of zeros of 04 and of sums of two of the dj ■ 
Zeros of a 4 : They are of the form fi = i(k + |)/z, k = 0, 1, 2, .... 
Zeros of 03 + a 4 ± : They are contained in the region r 34 ± : 

1 7T 

(Im/i) In 1 — - = Im £3,4 (/i) + Y(/i) — —Re /1 ± 7rRe \i. 

ImI 2 

Similarly the zeros of a 2 + a 4 ± are contained in T 2 4 ± : 

1 7T 

(Im /i) In - — r = Im Si^in) + — —Re /1 ± 7rRe /i. 
I A* I 2 



•7) 



The zeros of a\ + a 3 are contained in ri 3: 



1 7T 

(Im/i) In — = ImSi, 2 (/i) + -Re/t + F(/x). 
I A* I 2 



•9) 



The zeros of a± + a 2 are contained in T12' 



1 7T 

(Im/i) In — = Im5 3 All) + —Re u + F(/i). 

m 2 
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The zeros of a\ + a 4 ± are contained in r x 4 ± : 



11 7T 

(Im/i) In — = -(Im Si a + Im S 3A ) ± -Re/i + F(/i). 



.11) 



Put 



X = -Re/i + (Re/i)arg(-) - Im/i + ^ReO-(-) = -Re /i + Y~(/j). 
2 2 /i 2 



When Re /i > 0, a 4 + dominates over a 4 - and we shall only consider r 34 + = IY2, r 2j4 + = ri.3, 
r li4+ given by 

imS$ 4 + x, on r 3 4 + = ri 2 



(Im/i) In ji- = ^ Im5i, 2 + X on r 2>4 + = Ti^, 

IW I UlmS lt2 + ImS 3A ) + X onI\ )4 +. 



.12) 



Recall now from Section 5 that 0-(hj /i) in H8.6|) appears as a remainder in Stirling's formula, 
so that duO-Qi/ fi) — 0(h//j, 2 ), and hence X is uniformly Lipschitz for |/i| > h. Proposition lB.il 
can therefore be applied to get the approximate behavior of the T^. This will be exploited later. 

In the left half-plane, a 4 - dominates over a 4 + and we consider all the 5 curves r 3j 4- , r 2 ,4- , 
Ti,4-, r lj2 , and T li3 , given by: 



(Im /i) In 1 — - 
I A* I 



' Im 5*3 j4 — 27rRe /i + X, on T 3 4 - 
Im Si j2 — 27rRe /i + X, on T 2 4 - 
Im Sip + X on rj 3, 
Im 53,4 + X on Ti. 2l 

i(Im 5i i2 + Im 53,4) — 7rRe /i + Ion I\ 4- 



(8.13) 
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Again Proposition IB . 1| can be applied to give the approximate shape of Tjfc. Recall that we are in 
the case 1 with \p\ 3> h, so that (|8.2|) holds. 

1. Skeleton in the region Rep > 0. (We will implicitly use that |a 2 ||a 3 | = |a 1 ||a 4 +|.) The 
region |Re/Lt| < O(h) will require a special discussion. In the region Re p > Ch, we have \a£\ = 
e 2irRe/i/h| a4 _ | \a±\, so | a.4 1 ~ \a^\ and in this region we see from an earlier observation that 
the zeros of a 2 + 03 will not play any essential role. In this region we shall therefore use the 
appearing in (j8.12() and, as pointed out earlier, we are interested here in the part of each such Tj t k, 
where \aj\ = \dk\ dominates over the other \a v \. 

It follows from Proposition IB . 1 1 that the curves in (|8.12|l (as well as the ones in l|8.13|l ). are of 
the form 

Imp = lj,k(Ren), with | 7ii fc |, \i h k \ < 1. (8.14) 

Notice that every crossing point of two of the curves T 1 4+, Ti, 2,1^2 4+ is a l so a crossing point for 
all three. This follows from (|8.12|) or even more trivially from the observation that 2 of the 3 
equations |oi| = |flt4+|, |oi| = lag], | eta j = |&4+|, imply the third one. Also notice that if we draw 
the two curves Ti^ = 1^3.4+, Ti^ = r 2 .4+, then r lj4 + is between the two — see Figure 2 on page 35. 

For p > we have X(p) = hReO-(^) = hO{e~ 2 ' K ^ h ) by (JS^Jl and hence j jtk (Rep) is 
described as in Proposition IB . II with 



(ImS 3A (Rep), (j, fc) = (3,4+), (1,2), 
F = F j<k (Rep) = hO(e- 2 * Re ^ h )+ { ImSi, 2 (Re/i), U, k) = (2, 4+), (1,3), (8.15) 

[ |(ImS 1>2 +ImS 3 ,4)(Re/x), (j,k) = (1,4+). 

In the region Imp < min(7 2 ,4+ (Re //), 73 4 + (Re /1)) we have 1 0,4+ 1 > max(|ai|, |oa|, \as\, I04- 1), 
and if we restrict further to 

Ch 

Imp < min(7 2 4 +(Re/i),7 3 4 +(Re/z)) p, Re p > Ch, (8.16) 

In f-r 

with C»l, we see that 04+ is dominating in the sense that 

|o 4 +| > 2\ai +a 2 + a 3 + o 4 -|, (8.17) 

and hence G(p; h) has no zeros in that region. Similarly in the region Im p > max(7i j 2, 71,3) (Re p), 
we have |ai| > |aa|, 1 0.3 1 , |a 4 ± |, and if 

Ch 

Imp > max(7i 2,71 3)(Re p) + - — =-, Rep>0, (8.18) 

In -t-t 

with C»l, then a\ is dominating in the sense that 

|a x | >2|a 2 + a 3 + a 4 |, (8.19) 

and again G(/i; /i) has no zeros there. 

Now consider a point p £ T 3i + , where 734+ < 72,4+ , so that Im £3,4 < Im Si t 2, with Re ^> h. 
Going down (ie. decreasing Imp while keeping Rep constant) by a distance 3> /i/ In jiy , we reach 
the region (|8.16|l . where 04+ dominates. 

Case a. < ImSi :2 — ImS^ < C)(/i). Going up by a distance 3> ft/ln^r, we cross r 2:4 + = I\3 
and reach the region (|8.18|) . where a\ is dominating. 
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Case b. Im Si )2 — Im 5*3,4 > Ch for C>1. Going up by a distance ~ h/ In we reach the region, 
where 03 is dominating: 

|o 3 | >2|ai+a 2 + a 4 |, (8.20) 
and continuing to go up, 03 remains dominating until we reach a hj In y^y-neighborhood of IY3 = 
F 2 4+. After crossing that curve and going up by another amount Chi In jK, we reach the region, 

' IMI 

where a\ is dominating. 

Our discussion shows the following 

Proposition 8.1 We work in the region 18. 2|) and assume in addition that Re/i > 0. 
Iflmfj, < min(7 2i 4+,73 i 4+)(Re^), then 

|o 4 +| > max(|ai|, |a 2 |, |a 3 |, |a 4 -|), (8.21) 

// Im /i > max(7 12 , 7i.3)(Re then 

\ai\ >max(|a 2 |,|a 3 |,K±|)- (8.22) 

^/73,4+(Re/i)(= 7i l2 (Re/z)) < Im/j < j lt3 (Refi)(= 7 2) 4+(Re/i)), t/ien 

|a 3 | > max(|ai|, |a 2 |, |a 4 ±|). (8.23) 

^f72,4+(R-e^i)(= 7i,3(Re/x)) < Im/i < 7^4+ (Re /i) (= 7i, 2 (Re/z)), t/ien 

|a 2 | > max(|ai|, \a 3 \, |a 4 ±|). (8.24) 

If the distance from /i io (T± t 2 = r 3 4 +) U (T± t 2 = r 3 4 +) is > Ch/\njj^, with C 1, then in 
the respective cases (|8.22|) - (18.24(1 can be sharpened to the dominance in the sense explained above. 
In particular, G has no zeros in this region. If in addition, Re fi > Ch with C> 1, then we have 
the same conclusion in the case of (|8.21l) . 

In the region l|8.2[) . intersected with the right half-plane Re /1 > we define the skeleton to 
be the union of the curves Im\i — max(7i :2 , 7i i 3)(Re fi) and Im pL = min(7 2 4 + , 7 3 4 +)(Re /i). The 
proposition shows that the zeros of G in the region under consideration are contained in the union 
of all discs D(/i, Ch/ In y^y) with /i in the skeleton just defined, and the set of all pi below the 
skeleton, with < Re \i < Ch, for C>1. 

2. Skeleton in the region Re/i < 0. Again the region |Re/Lt| < O(h) will require a separate 
discussion so we restrict the attention to Re^ < —Ch and we will use |a 4 -| = e~ 27Ttl / h \a 4 + | ^> |a 4+ |, 
so that \ai\ ~ |a 4 - |. We therefore concentrate the attention to the curves in (|8.13|) . As before we 
notice that every crossing point of two of the three curves I\ 2 , T 2 4 - , T 1 4 - is a crossing point of 
all three. The same holds for IY 3 , T 3 4 - , T 1A - . 

Now use that is small and hence that Im Si 2, Im S^^ and their derivatives with respect 

to Re /j, are small. We can therefore consider the two curves: 

A: — 27rRe pi = Im Si j2 — Im ft, 4, B : — 27rRe /i = Im 5*3,4 — Im Si, 2 . 



3G 



Im /i 




Re fj, 



Figure 4: The union of the solid curves in the figure gives a schematic representation of the skeleton 
S' in the right half-plane intersected with the region 1)8. 2[1 . Proposition 8.1 shows that the zeros 
of G in this region are inside the union of the thickened skeleton, obtained by placing a disc of 
radius Chj |ln|/i|| around each point fi E S' , and the set of all /i with < Re/i < Ch below S'. 
Proposition 8.3 gives a more precise description of the location of the zeros of G with |Re fi\ < O(h). 
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They are of the form 

— Re M = 7a (Im /x), — Re fi = 73 (Im /x), 
where 7,4, 7s are small with small derivatives and satisfy: 

|7A(Re/x)| - |7 B (Re/i)|, 7 a7b < 0. (8.25) 

The curve I\ 4 - will play a central role. It crosses A,B at the points MAjMb (unless these 
points arc hidden in the forbidden region), and we have 

(Re //a) (Re < 0, \Refi A \ ~ |Re/i B |. (8.26) 

We notice that ma is the unique crossing point for I1.3, r 3 4 - , r 14 - while /is is the unique crossing 
point for 1^2, r 2;4 - , r 1;4 - . More precisely, 71,3(0 — 71,4- (t), 71,4- (i) — 73,4- (0 vanish precisely 
for t = Re \xa and have the same sign as t — Re ha- Similarly 71,2 (i) — 7i.4- 7i,4- (i) — 72,4- (t) 
vanish precisely for t ~ Re mb and have the same sign as t — Re /is- We also notice that if ji 
belongs to one of the three curves I\ 3 , Ti 4- , T 3 4- , then the distance from /i to any of the two 
other curves among these three is > C _1 |Re/i — Re ma | / In y^y . The same observation holds for 
Ti.2, f\ 4 - , r 2: 4- with mb instead of ma- 

For m < 0, we have X - vrRe m = /iReO_(^) = /?.C(e^ 27r l^l/ /l ) by and hence 7l 4 - is 

described as in Proposition IB . II with 

F = F 1A -{Rey) = ^(ImS h2 + lmS 3A )(Rey) + hO{e- 27rlRe ^ /h ). (8.27) 

Now assume to fix the ideas that Re ma < (the case Re/iB < can be treated similarly). 
Considering the three curves IY3, r 14 - , r 34 - , we see that for Re m < 0: 

If Im/i < min(7 lj4 - , 73 4- )(Re /i), then 

1 04- I > Kl, l a 2|, \d3\, |fl4+|- 

(In this case, we also have Im y < 72,4- (Re /i)-) 
If Im [i > max(7 lj4 - , 7i j3 )(Re m), then 

|ai| > M, \a 3 \, \a A ±\. 

(In this case, we also have Imy > 7i j2 (Re/x)-) 
If 73,4- (Re //) < Iti m < 71,3 (Re m), then 

|as| > I ai|, |a 2 |, |a 4 ±|. 

This covers all possible cases with Re/i < 0. Notice that the last case can appear only when 
Re ma < Re M < 0. When Re Ms < 0, we get the analogous discussion after a permutation of the 
indices 2 and 3: 

If Im/i < min(7i 4- , 72 4- ) (Re m) > then 

1 04- 1 > |oi|, |a 2 |, |o 3 |, |o 4 +|. 
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If Imfi > max(7 14 - , 7i i2 )(Re /x), then 

Kl > I02I, \aa\, |o 4 ±|. 



If 72,4- (Re /t) < Im /1 < 71,2 (Re /i), then 

M > Kl> l«3 U l«4± I 



Moreover, if in addition all the inequalities for Im/t are valid with an extra margin C/i/ln-r^r, 
C> 1, then in the various cases, we have dominance of a%, 02, 03 respectively, in the sense explained 
before. For the dominance of 04- , we also need the assumption that Re/x < —Ch with C>1, 

Exponential localization to the skeleton. Recall that we are still working in the region l|8.2[l . In this 
region we define the skeleton to be: 

S = S"ur 4 , (8.28) 

where we define S' to be the union of the following two sets in respectively the left and the right 
half-planes: 

- In the closed left half-plane (intersected with (|8.2[l ). assume Re HA < to fix the ideas, then this 
part of S' is given by all points of the form Im/x = 7 14 -(Re/i) with Re /t < Re ha, all points of 
the form Im /i — 73,4- (Re /t) or of the form Im fi = 71,3 (Re /t), with Re fiA < Re fi < 0. 

- In the closed right half-plane the corresponding part of S' is defined to be the union of the two 
curves: Im/i = max(7i ! 3, 71,2) (Re/x) and Im/i = min(72 ) 4+, 73,4+) (Re/i). 

T4 is defined to be the part of the imaginary axis given by 

Re/i = 0, Ch < Im/i < min(72 j 4+,7 3; 4+)(0), 

where C is the same constant as in l|8.2|l . Notice that this part may be empty. The earlier discussion 
shows that we have 

Proposition 8.2 The zeros of G in the domain l|8.2|l are contained in the set 

|J D I n, -^j- J [J{fx below S'; |Re /Lt| < Ch}. (8.29) 

Remark. The localization result of Proposition 8.2 improves if we are far from the branch points 
of the skeleton. Thus for instance, if ju is below S' and dist (/t, S') > Ch/ ln t-t , then 



max 



|a 4± | >e^ dist ^ s ')(| ai | + | a2 | + | a3 |), 



so the zeros of G are exponentially small perturbations of those of 04 — a 4 + + a 4 - : In this region 
there is a bijection b between the zeros of a 4 and those of G, with 

|&(/i) - /i| < C(/i)cxp [-(C^r 1 ln(-^)dist (/t, S')]. 

I A* I 
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Im fi 




Re /i 



Figure 5: The union of the solid curves in the figure gives a schematic representation of the skeleton 
S' in the left half-plane intersected with the region i|8.2|) . Proposition 8.2 shows that the zeros of 
G are inside the union of the thickened skeleton, obtained by placing a disc of radius Ch/ |ln \pb\\ 
around each point /.i £ S', and the set of all \i with |Re^| < Ch below S' . 
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Similarly let fi be a point of S' in the right half-plane, say fi € I\ 3 with ai,a 3 dominating 
above and below this point respectively. If 

C 2 ft 

dist(/xo,r 34+ =Ti 2 ) > j— , 

In T-^-r 

then in D(/j, , Cft/(ln(l/|//|o|))), we have 

max(|ai|,|a 3 |) > e (chrl(ln ^ )dist(w ' ri ' 2) max(|a 2 |, |a 4 |), 

and we conclude that the zeros of G are exponentially close to those of a\ + a 3 . Essentially the same 
statement holds when /i belongs to the lower part of S', but here the size of Re \x also matters, so 
near a point no & ^3.4+ H S' we get a bijection & between the zeros of a 3 + a 4 + and those of G with 

6(/i) - /i = 0(l)(-V)(e-^ dist( ^ ri ' 3) + e" 3 ^), 

I Mo I 

and we have to assume both that dist (/io, ^1,3) 3> ft/ In t^t and that |Re /i| » ft. 
The analogous statements hold in the left half-plane. 

More refined analysis in the region |Re/i| = 0(h). The study of the upper part (Ti, 3 or Ti^) of 
the skeleton is unchanged in this region, while the lower part requires more attention in view of 
the fact that |a 4 | may be considerably smaller than max(|a 4 +|, |a 4 -|) when we are close to a zero 
of a 4 . In order to fix the ideas we assume that 7 3 . 4 +(0) < 72,4+ (0). 

After multiplication of a 4 , a 4 ± , a 3 by the same exponential factor, we arrive at 

ft 4 = a 4 + +a 4 - = 2cosh^, a 4 ± = e^, a 3 = e ^ 4 '^' h) , (8.30) 

h 

and we shall drop the tildes in the following discussion. Here 

<HM; ft) = S 3A (j*) ihO-(-) + 11 In % + ^ - n, (8.31) 

with 

- lm(f>(fr, ft) = -lmS 3 4 (/z) + Re — - ftRe £>_(-) + Im/iln-i- - Re^arg(-) + Imfj,. (8.32) 

2 /i 1 

We have 

& m „(-lm0) = In ^ + 0(1), 

0Re„(-Im0) - 0(1), 

V«(-W) = 0(M 1_H ). M>2, 

and recall that we work in the region \fj,\ > Ch. Notice that ftV M ln|a 3 | = V M (— Im 4>). Similarly, 
we look at 

sinh^ / 1 

^^'^^"(^W' (8 ' 33) 
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for a suitable branch of the square root. Also /i9 fr lna 4 = 0, so this relation gives a bound for 
/iV M lna4 and its real part; /iV At ln|a 4 |. We have the general estimate 



so we get 
Assuming 

we conclude that 
and consequently, 



|coshz| > idist(z,cosh" 1 (0))e |Rez| , for \Rez\ < Const., (8.34) 

7r|Re fi\ 

hd * ta - - ' sg " (Re "> + o(i) d is ,(4 : .^-'(o)) ^ (8 ' 35) 

dist (fi, - cosh _1 (0)) > — with C> 1, (8.36) 
7r In A- 

m 

|W At ln|a 4 || </i|V M ln|a 3 ||, (8.37) 



^i mM ln^4 = (l + (l))ln^+0(l), hd Re ^lnp\=0(l)+o(l)ln-^, (8.38) 
|a 4 | | A*| | <3-4 1 IMI 

where 0(1) denotes terms that are uniformly bounded and o(l) denotes terms that tend to 0, when 

dist U, ^cosh _1 {0}) 
— — > oo 

For each zero fij of cosh ^ , we introduce the diamond shaped neighborhood 

Ch 

Dj = {/i; |Re yu| + — Imud < (8.39) 

In A 
M 

with C large enough so that the preceding estimates apply away from the union of all the Dj. 
Define r 3i 4 to be the set of points with |a3 1 / 1 0.4 1 = 1 away from the union of all the Dj, with Dj 
added, if Dj has the property that the distance from this diamond to the points just defined, is 
zero. Dj is unique if it exists, since the other points of r 3i 4 form a curve \m\i — 7 3i 4(Re 11), with 
1 73 4 1 <C 1. From the above estimates we get 

In In A O(h) 
l3A± (R ef i)-0{h)—-ft < 73,4^) < 73i4 ±(ReM) + r^r-, ±Re M > 0. (8.40) 

m t— t hi A 

M I Ml 

(If T 3j 4 stays away from an /i/C-neighborhood of the zeros of cosh then the agreement is better: 

73,4 ( Re M) =73.4±(Re/i) + r^p, ±Re^ > 0.) 

In A 

I Ml 

In fact, we can get an even more precise estimate for the distance between T 3i 4 and r 3 4 ±: Let 
/i € r 3j4 , and put 

d(fM)) =max(- j— ,dist (fi , oJ^O))). 

In At 

I mo I 
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Then 



d(jM>) < a4(Mo)l < ^ 



C7i |a 4 ±(/i )| 

away from the diamonds. We therefore get the following estimate for the vertical distance from /xo 

to r M ±: 

|73, 4 (Re^ ) -7 3i4 ±(Re/i )| < In - J7 - T , (8.41) 



ln^r 

I Mo I 



d(/i ) ' 



(assuming for simplicity that rf(/io) < ft/2). This is a refinement of the lower bound in 1)8. 40[) . and 
the argument also gives the upper bound there. 

We reach the following conclusion about the location of the zeros in the region |Re/x| < Ch: 

Proposition 8.3 - Above S' and at distance > Ch/ In from S' , ai(fi) is dominating. 

- a4 is dominating if [i is below S' , at distance > Ch/ In yiy from a 4 ~ 1 (0) and at distance > 
Ch^{^) from S' . 

- In between (for instance below but above r 34 ± ), a 3 (or 02) is dominating if the distance to 
the skeleton is > Ch/ hii—r. 



Improvement in the region Re fi» h. Let us recall that 

aia 4 + = a,2a 3 . 

Therefore, 

a\ + a.2 + a 3 + ft 4 + = a 4 + ^1 
The zeros of 1 + (fl2/a 4 +) are situated on T 2i4 + and are given by the explicit quantization condition 

(8.44) 



Q2 
a 4 + 



1 



a-3 
a 4 + 



(8.42) 
(8.43) 



Ulnu - u + — + Sx 2 + ihO(-) = 2ith(k + -), fee Z. 
4 /1 2 



The distance between successive zeros is ~ h/ln-^. If [iq is such a zero, then in a disc D(fj,o,r) 
with r <C ft / In -i-J-r , we have 

' I Mo I ' 

a 2 



a 4+ 

Away from the union of all such discs, we have 

a 4 + 



Im-mo| 



ImoI 



> 



0(1) 



Similarly, the zeros of 1 + (a 3 /a 4 +) are situated on the curve T 3 4 + and given by the quantization 
condition 

uln/i- fi + + S 34 + ihO-(-) = 2irh(k + -), fee Z, 
4 [i 2 

and the other statements about 1 + (a2/a 4 +) carry over to 1 + (a 3 /a 4 +). 
Now consider 



.45) 



G([i; h) = a 4 4 
We get 



1 



«2 

a 4 + 



1 



Q3 
a 4 + 



a 4 - 
a 4 + 



= a 4 + 



1 



_&2_ 

a 4 + 



1 



Q3 
a 4 + 



(8.46) 
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Proposition 8.4 In the region Re fi 3> h, there is a bijection b from the union of the zeros of 
1 + ^2 / &4+ o,nd of 1 + as /a 4 + to the zeros of G with 



(8.47) 



So in the region Re /i» h, and modulo an exponentially small error, we can identify the zeros 
of G with the union of the zeros of 1 + 02/04+ and of 1 + 03/04+. 
This finishes the analysis of the skeleton in the first case 1|8.2|) . 

Case 2: Assume that 



arg/J + - 



< 7T- i /i< |mI « 1. 



(8.48) 



In this case from (|5.8|l and l|5.9|l . respectively, we recall that 
a2,3 = 2 cosh (~r~) ex P 



,h . i , . , . 77/1 „ 

(j, h 4 



7T/i 

2ft 



and 



Using |JOJ, we get 



01,4 = exp 



-0 + (-)-^(MlnM-/,+ T ) + ^ 



F(/z; h) = e 2h G([m; h), G(p; h) — a\ + 02 + 03 + a 4 , 



(8.49) 



where aj, j = 1, 2, 3, 4 are the same as in case 1, but now with a partially different representation: 



Oi± 



aj+ + a 1 - , 

e i(Si,2+S3,4)+0 + (£) + £(Mln(v)-M+^)±¥ 



(8.50) 



a 2 = e iSi,2+{ 



«.3 

04 



e <> 



53,4+4 



-O+(£)-jr0*M«M)-M+ 3 £) 



Again we consider /i times the real parts of the different exponents of the aj, j = 1 , 2, 3, 4 
ri± = — Im Si 2 — I™ S3 4 + (Iw /i) In - — - — Y (/1) ± 7rRe fj,, 

ImI 

r 2 = -Im5 L2 + -Re ^, 

7T 

r 3 = -ImS^ + -Re/J 

r 4 = -(l777 /j)ln-i- + ?(/j), 

ImI 



where 



F(/j) = (Re /j)arg (i/j) — Im/j — ftReO + ( — ). 



(8.51) 
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By a symmetry argument, we shall now see that (^±^2,7*3^4) plays the same role in the 
present case 2 as (r 4 ± , r 3l r%, r\) in case 1, provided that we perform the following transformations 
on (r 1 ±,r 2 ,r 3 ,r 4 ): 

1) Add Im (5i,2 + $3,4) to each of the five terms. 

2) Replace /x by ~p. Then we get rj(~p) = (Im (Si,2 + £3,4) + rj)(/j,): 

r 1± (p) = -(Im/l)ln T 4r + (Re u)arg(?) -Imp + /iRe£> + (-) ±7rRep, (8.52) 
I A*l * M 

r 2 (/i) = ImS 3 , 4 + -Rept, 

_ 7T 

r 3 (fJ-) = lmS 1:2 + -Re/j,, 

r^Cp) = Im S\ 2 + Im S3 4 + (Imp) In — — — (Re p)arg ( — ) + Imp — hRe 0+( — ). 

This is analogous with (r4± , r3, T2, ri) in case 1 except for the fact that ReO + (^) here corre- 
sponds to ReO_(^) in case 1. 

Remark. Using Q5.5JI . it is easy to check that 

?(/*) - Y{p) = ±7rRe fx + hO{e- 2 ^ Ke ^ /h ), (8.53) 

when 

|arg(± M )|<|-i | M |>ft. (8.54) 

It follows that 

= r J ±( f i) + hO(e- 2 ^ Retl ^ h ), j = 1,4, (8.55) 

when ll satisfies Q8.54p . and hence if fx belongs to the skeleton S' defined according to Case 1, the 
distance from ll to the corresponding skeleton S' defined according to Case 2 is 



h 



In X 

m 



-27r|Re/x|//i 



We end this section by some general considerations that will be useful in Sectional We see 
from H8.15J) that the spectrum will have a genuinely 2-dimensional structure if 

|ImS 3 ,4(0) -ImSi, a (0)| > h, (8.56) 
or if Im 63,4(0) and Im £1,2(0) have the same sign and 

min(|Im^ 3 4 (0)|, |lm^ 12 (0)|) > h\n\. (8.57) 

h 

In the latter case, we even have some eigenvalues on the imaginary /x-axis, related to 1-dimensional 
barrier top resonances. It is therefore important to have a sufficiently invariant and direct descrip- 
tion of Im 53,4(0), ImSi :2 (0). 
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The final definition of in the beginning of section 6 is simply that we start with the null- 
solution /4 of Q near «4 and extend it along the exterior part of -Ko,o until we reach a neighborhood 
of 0:3, where we get exp (^ 53,4)73. (Here we neglected the real Floquet parameter 0\, since we are 
only interested in the imaginary part of 53,4). The definition of is similar. 

Now take n = (cf. (|6.14|1 and l|6.15[l ) and represent the operator Q as acting in H^ C (Q), where 
<E> is strictly plurisubharmonic, with A$ ~ a neighborhood of Kq.q in T*R. From the construction 
of ej, fj, we see that fj is near ctj a normalized null-solution of Q in H$(Cl), where $ — $0 is 
small and $0 is defined in a sufficiently large neighborhood of the projection of the branching point. 
Here if kt is the canonical transformation associated to some standard FBI-Bargmann transform, 
then A$ = k t (R 2 ), A$ = n T o ^(R 2 ), with £/ as in (|6~Ti|) . Since A$ = {£ = 7^-} is naturally 
identified with T*R, where p = x£, so (since \i — 0), we know that the null set of Q intersects 
A$ along two crossing "real" curves that we we can identify with "the interior part" of K . 

Extend $ to be defined in Q with $ - $ small. If 5 3j4 ~ S° 4 + /i5^ 4 + then 

- Im Sg )4 = / (-Im (£ • da;) - d$ ), (8.58) 
•^73,4 

where 73,4 now (cf 1)6.24(0 is a real curve from to 03 in Qq 1 (z(0)) close to the exterior part of 
the "left loop" of i^o.o- Here we let z(0) denote the z-value in (|6.15|) corresponding to /1 = 0. Let 
this left loop be denoted by 7 1 and let us consider it (after slight deformation) as a closed curve in 
Qq 1 (z(0)) joining the critical point of Qo to itself staying close to the left loop of K . Here, we 
may assume that the interior part of 7 1 (joining 0:3 to 0:4) is contained in A$ , so — Im (£ • dx) = e?$o 
there. Hence (|8.58|l becomes 

Im5g 4 = / (Im (£ • dx) + d$ ) = / (Im£-dx). (8.59) 

Here Im (£ • da:) in the last integral can be replaced by Im (£ • dx) + d$, which by Stokes' formula 
can be further replaced by any other real 1-form lo with du = Ima, cj| A = 0. This means that we 
can reinterpret the last integral in H8.59|l as the corresponding one along the corresponding closed 
curve in the complexification of R 2 . 

To simplify things further, recall that e, h 2 /e are small perturbative parameters for Qq and that 
Qq 1 (z(0)) is real when e = h 2 /e = 0. In general, if q = q s depends smoothly on a real parameter 
s, if E is not a critical value and 7 = 7(s, E) is a simple closed curve in q^ 1 (E), then for E fixed; 

§-J tdx = -£ E ' S) |f (*(*),£(*))<**, (8.60) 

where [0, T{E, s)] 9 t 1— > exp i-ff 9s (p(0)) is a natural parametrization. This can be applied to the 
case q s = q — E(s), so if E also depends on s, we get 

s/,^-r^-^«')* 

In this form, we can treat a loop like 71 (s) C q,r 1 (-^'( s ))j starting and ending at the critical points 
p c {s) of q s , parametrized by ] — 00, +oo[i— * exp (tH qs )(p(0)), provided that we take E(s) equal to 
the critical value q s (p c (s)): 

§-J Jdx = £™ (J-(q s (p c ( s ))) - (^q a )(x(t),t(t))j (8-61) 
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This can be proved by a limiting procedure, approaching 71 (s) by closed curves at non-critical 
levels. 

Taking the imaginary parts, this means that we have a fairly simple way of computing Im S\ t 2, 
Im 53,4 perturbatively. From this computation, we see that it is of interest to compute the e 2 
contribution to the averaged principal symbol. This computation was carried out in Section 2 
of under the assumption that (q) = and it works essentially the same way without that 
assumption: We start with the principal symbol 

p f =p + ieq + f 2 r + C(e 3 ). (8.62) 

The function 

1 pT(E) / T(E)\ 

Go = j [t - -^i J q o cxp (tH p )dt, on P ~ l {E) (8.63) 

introduced in Proposition 12. II satisfies 

H P G = q- (q). (8.64) 

Put G = Go + ieGi + 0(e 2 ), where G\ remains to be determined. As in p|] we get at a general 
point cxp (ieH G )(p) € A eG , (p e T*M): 

Pe| Ae ^ p e (exp (ieH G )(p)) = ^ €l [f Pc(P) 



= p + ie(q) + e 2 (r + H p Gi - H Ga {\{q + (?)))) + <D(e 3 ), 

where we used that H Go p = —H Go (q — (q)). 
Letting G\ solve 



H p G x = H Ga {\{q + (q))) - (H Go (\(q+(q)))) - (r - (r», 



we get with G = Go + ieGi, 

p e (exp (ieH G )(p)) =p + ie(q) + e 2 (s) + 0(e 3 ). (8.65) 
Now assume for simplicity that T(E) — T is constant. Then 



(s) = (r) 
= (r) 



1 f T 

— {G ,q+{q)}°exp(tH p )dt 
1 f T f T T 

T 2 J J ^{ qoex v( t + s ^ H P'( q+ w) o exp {tH p )}dtds 



r-T i-T 

2T 2 . 

1 /" T , T 



( s ~ -)({? oex P (g>})ds. 



w 2T./n 2 

Here, we notice that 

({goexp ( Sj ff p ),< g )}) = {<goexp( Sj H p ),( g )} = {< g ),(g)}=0, 
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so finally: 

( S )=(r)-iy ( s ~^)({qoe X p( s H p ),q})d s . (8.66) 

The formulas (|8.65|) . (|8.66|) will be used in Section ^| together with the following remark: If 
we put 

Cor (q 1 ,q 2 ;s) = ({qx o exp (sH p ), q 2 }) , (8.67) 
then a simple computation shows that 

Cor (q 1 ,q 2 ;s) = -Cor (q 2 , g x ; -s). (8.68) 

If we put 

1 f T , T, 



C( qi ,q 2 ) = -^ (s--)Cor(qi,q 2] s)ds, (8.69) 



then combining (|8.68|) and the T periodicity of Cor (qi, q 2 ; s) with the change of variables T/2 — s = 
s-T/2, we get 

C(q 1 ,q 2 ) = C(q 2 ,q 1 ). (8.70) 

9 Skeleton for \/i\ < 0{h) 

In this section we shall consider the case < 0(l)h. In doing so, we will use (|4.13|) more directly. 
Case 1. We will first work in a region {/i G C; |^| < rh} U {/i € C \ {0}; |arg/x — || < ir — 1/C}, 
where < r < 1/2, C > 0. (The corresponding region with |arg/z + f | < ft — 1/C, can be treated 
with a symmetry argument as in the end of Section [SJ and this argument will be given later.) It 
follows from (|4. 1 3|> that here 02,3 ^ and lnT(i — i^) is well-defined, while T(i + may have 

zeros. Consequently we use the reflection identity, to get 



V 27r u- 

aiA = WW) 



2ir h ' 



Now using H8.1(l we get 



In F( 2 'ft' 
+ e v^F 

e^G(fi;h), 



e '" F ' ' +' £ - t + % +1" 2 cosh ay 



where 



G(/i;ft) = e^+^-^l^-lnl^+if (9.1) 



+ e h h ^ v^f 4 2 cosh — . 

h 

ai + a 2 + a 3 + 04, a 4 = a 4 + + a 4 - , 
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where the terms are the same as in (|8.3|l , although we shall now use different asymptotic approxi- 
mations. 

Again we introduce h times the real parts of the different exponents: 



n = 


— ImSxfl - 


Im S3, 4 4 (Im /Li) In — 


— hRe In 


r<i = 


— Im Si, 2 4 


7T 

-Re/i, 




rz = 


— Ito5 3 , 4 4 


7T 

-Re/i, 




r 4 ± = 


— (Im /i) In 


) + hRe ln T ^Zl^ 
h 


± 7rRe /i. 



r(|-^) 



2tt 



(9.2) 



As before, we have 



T2 4 r 3 = ri + r 4 + . 



(9.3) 



Again, we define the different curves Tj.k by \aj\ = |ofc|, for j k £ {1,2,3,44} with the 
exception of (j, fc) = (2,3) and (j,k) = (4 + ,4~). (The segment 14 is now defined to be the 
segment of the positive imaginary axis joining to the lower part of S' , provided that this lower 
part is not hidden in the forbidden region, in which case we let 14 be empty). More explicitly, we 
get: 

Im S3. 4 ± 7rRe [i — ^Re /i 4 A on r 3j4 ± 
■ Im Si. 2 ± 7rRe [i — ^Re /i 4 A on T 2 ,4± 
(Im/i)ln- = < ImSi, 2 4 f Re/i 4 A on 14,3, (9.4) 
/! 1 ImS 3 , 4 4 §Re/i4 A on ri, 2 , 

i (Im Si, 2 4 Im S 3 , 4 ) ± S£ 4 A, on r 1)4 ± , 



with 



A = hRe In 



r(|-^) 



/2tt 



(9.5) 



The function A now differs from that of Section[S]by a term — -|Re /i. The definition of 14. 3 
14,2 = r 3j4 + coincides with that in Section [S] in the overlap region. 

We shall also define a set Tj A for j — 1,2,3 as in the preceding section. To do so, we check 
that 

1 



h 



Vcosh^ 



cosh^ 



< In- 



if 



In this region, we also have 



dist (fi, - cosh 1 (0))>-^ T . 
7r In i 



0(1) In 

In the region l|9.7|) we can define Tja by \a 



< 



a 4 



a4± 



<0(1). 



(9.6) 
(9.7) 

(9.8) 



1 0,4 j , and see that we get a curve 



Im /i = 7^4 (Re /i), 



(9.9) 



4!) 



with |t£ 4 1 <C 1. Using (|9.8|) . we also see that if we represent Tj 4 ± by Im = 7^,4* (Re fi), then 

(9.10) 



( 7j . 4 - 7j . 4± )(Re M ) = 0(l)/ l (^)(i). 



Actually the upper bound can here be improved to 0(ft)/ln^-. In analogy with Section [H] we 
define a diamond shaped neighborhood of each zero [ij of 04 by 



Dj = < (j,; |Re fi\ + \lm fi — Ira fij \ < — r 1 , 
I ln 7I J 



(9.11) 



The previously defined Tj.4 can hit at most one of the D v and if that happens, we add that diamond 
to the set 1^4. 

Now define the skeleton as before: 

s = s'ur 4 , (9.12) 

and as before we can describe the regions of dominance: 

04 dominates at distance ^> Mt^Xtt) below infj = i : 2,3 7j,4± , for ±Re/j, > 0, intersected with 
the complement of the union of the diamonds. 

The other cij dominate according to the earlier rules in their respective regions at a distance 
3> h/hxj- from the skeleton. 

Case 2. We now consider the case when /i belongs to the set 

1 



{/iGC; < rh} U {fi S C \ {0}; 
where < r < \, C > 0. From 14. 131 . we get 



arg fi 



(9.13) 



01,3 = -02,4 = e 



ai,4 



where 01,4 is non- vanishing, while 



'2tt 



02,3 



r(i-<«) 



ft'h g 2h + 4 



may have zeros, so we use the reflection identity to write 



02,3 = 



r (2 + % h) Si+is 



2tt 



ft* ft e iir+" 2 cosh ^— . 



We then use (|8.1|1 to get 
with 

ax± = exp 



F(»; h) = e^G((i;h), 



G(/x; ft) = a\ + 0,2 + a 3 + a 4 , ai = 01+ + ai 



x (5 1>2 + 5 3l4 ) + In! 



17T 7T/i 



if lnft+^i^f 
ft 4ft 
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ai = exp 



2^ 



03 = exp 



h S3A+ 2h, 



y ' h 4 



04 = exp 



Again, we introduce h times the real parts of the different exponents: 



r x ± = -Im 5i, 2 - Im 5 3)4 + (Im (i) In - + /iRe ln( ^- ) ± TrRe /i, 



7T 7T 

r 2 = — ImSxp + ^Re/i, r 3 = -ImS 3i4 + — Re/z, 

r 4 = -(Imu)ln- - hRe ln(- 2 h ^ 
h 



2ir 

We shall now make the same symmetry transformations as in Section 8, to see that the functions 
r i± * r 2, r 3, and ^4 play the same role as r 4 ± , r 3 , r 2 , and r*i respectively, in the previously considered 
case: 

1) Add Im (Si,2 + 53,4) to each of the rj. 

2) Consider the rj as functions of ~p. Then we get 



1 F(l-iH) 
r 1± (/z) = -(Irrc/l) In - + hRe ln( 2 _ fc ) ± TrRe/Z, 

« V 27T 



r 2 (//) = 177163,4 + 2^ e Mi 7 '3 = ItiSi,2 + ^Re^, 

_ 1 r(i-zS) 

r4 = Im 5i t 2 + Iwi 63,4 + (Im /I) In — — /iRe ln( — 2 - ). 

" v 27T 

Thus apart from a change of sign in ImSi^, ImS^, we see that (r"i± , F2, F3, r^) has the same 
properties as (r4±,7*3,r2,r'i) in the previously considered case. 
Remark. In the overlap region 



/)((). ,7,1 iJ ^ w | M | < C7i, |argM| < | - ^, or |arg(-/i)| < - - - 



where both cases apply, we notice that trivially rj = r 3 ± + 0(h) for ±Re fi > 0, j = 1,4, (and 
these estimates improve by (|8.55|l when \n\/h increases). As in the remark at the end of Section 
[SJthe distance between the two skeletons, defined according to Case 1 and according to Case 2, is 
therefore 0(h/ln^). 



10 Eigenvalue counting 

In each of the cases 1 and 2 of Sections |S| 03 we defined a skeleton 5 consisting of a horizontal 
part S' possibly with a vertical part (r 4 in Case 1 and Ti in Case 2) added. We notice that the 
definitions in the two sections agree for each of the cases 1 and 2 in the overlap regions for the two 



51 



sections, and we saw in the remarks at the end of the sections, that if we compare the skeletons 
for the two cases in the overlap region 



{|Re//| > -|Im/i|}U£>(O,r/0, 



then the distance between the corresponding skeletons is 



0( 



h 



ln(//)h 



e 



Now define the body by widening the skeleton: 




pes' \ 



In 



Ch 



1 



) 



UB„U5 e . 



(10.1) 



Here B Vl B e may be empty and will now be defined. They are non-empty if S' stays entirely in 
the admissible regions for one of the cases, and in order to fix the ideas, we assume that this is 
Case 1, and S' does not intersect the negative imaginary half axis. If so, we have a non-empty 
segment 1?4 in the imaginary axis, joining to the closest imaginary point of S' . Recall that we 
have defined the diamonds Dj around the zeros of 04 in I^, by (|8.39I9.11|I . We define B v to be 
the union of (in Case 1, and T\ in Case 2) and the corresponding diamonds. (In Case 2 we do 
the corresponding definition with "4" replaced by "1" and Ti is then the segment in the negative 
imaginary axis, joining to the closest part of S'. ) B e is non-empty precisely when 1?4 or Ti is. 
In Case 1, it is defined to be the set of points \i below S' at distance at most 



from S' with C > sufficiently large and with |Re/i| < h. Here the upper bound h in the last 
estimate, may be replaced by h/Co for any fixed Cq > 0, and we could decrease B e further by a 
more detailed discussion. In Case 2 we have the analogous definition. 

We next define what we mean by an admissible curve. It should be a piecewise C^-curve 
7 : [a, b] — ► C without self-intersections, parametrized by arc-length. It is tacitly assumed that we 
consider a family of such curves, which is uniformly bounded in the sense that we have uniform 
bounds on the number of jump discontinuities of 7, the continuity of 7 between the discontinuities, 
and on the length 6 — a. It is also required that 7(4) may belong to B only for t £ Ij, j = 1, 2, M, 
where Ij are disjoint intervals of length < Ch/ In j-^, for some /1 = j(t) € Ij, if j{Ij) H B e = 

and of length < Ch^^-(j-k-) otherwise. We also assume that we have a uniform bound on the 



number M of such intervals. 

Assume for simplicity that a, b $ U/j and let us partition [a, b] into intervals in increasing order: 



For each J^, let a v (u\, be the corresponding dominant term along 7(Jfc). For simplicity we shall 
assume that the image of 7 is entirely contained in the admissible region for one of the cases 1 or 
2, so that v(k) is either in {1,2,3,4*} (Case 1), or in {1*, 2,3,4} (Case 2). Let 




[a, b] = Jo U h U Ji U I 2 U ... U I M U J m . 



(10.2) 



Mfc,e = Mfc+M = 7(*fc+i) for some tk+i ^ 4+1, k = Q, ...M - 1 
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and put 

Mo, s = 7(a), MM,e = 
Then we have, with aj = e^^ h (cf. (|E3Jl, CUD): 

Theorem 10.1 Let 7 6e an admissible curve as above. Then 



1 f G' 1 



M-l ^ 

X! (<t>v(k)(Hk,e) - <t>v(k+i){Vk+i,s)) ~ <MMo, s )) + 0(1) + 0(maxlnln — — )) 

fc=0 



where the maximum is taken over all k with "/{Ik) H i? e ^ 0, so i/7 never meets B e , we only have 
the remainder 0(1). 

Proof: Notice that the first term of the right hand side of l|lU.3fl can also be written 

1 M 

Re -— 2_,{<t>u{k) (jtk,e) - <l>»{k)(Hk,s))- (10.4) 

k=a 

Consider an interval Jk- If we first assume v{k) ^ 4 ± (if we are in Case 1), then for t € Jk- 



G(l(t)) = 6, W (7W)a,( fc )(7W), \bu{k)(l(t)) - 1| < \. (10.5) 
Let jufc.s and JLk,e be the start and the end points of 7^ , so that 

MM = MM+0(^7^(7— V))> (!0-6) 
In (Mfe,s)h 

~ „ In In 1 

Mfc.e = Lik,e+0(h- — (- — )), 

m {Hk,e)h 

with the improving to if the corresponding neighboring interval I,, does not meet B e . Using 
<|1U.5[1 . we see that 

IPG' 

2ni -g-^M = °( 1 ) +^(fc)(Mfc,e) - ^(fc)(^fc, s )- (10-7) 

In the case v{k) — 4 ± , we know that a 4 is dominating along 7^ and i|lU.7[l holds with 4> u ^k) 
replaced by h In 04. Now we also know that along 7| j , we have 

a 4(7(*)) = c(7W)Oi/()t)(T(*)). 

with 

l/(ClnJ_)<|c(7(t))|<C, 

and with arg 0(7(4)) of bounded variation. It follows that the real part of the equation l|l().7(l still 
holds in this case. 
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We next estimate the integral along 71 , and let us consider the worst case, when j(Ik)P\B e ^ 0. 

' J k 

Let 

,„ . , In In . 1 , . 
r = 0(h— — )) 
m (fi)h 

be such that 

7(/ fe ) C D(Mfc-i,e, £)• 

On this disc, we write 

G(pi; h) = a„( A ._ 1 )(/x)& k (ji), (10.8) 

where 6fe is holomorphic, and 

C > |& fc (Mfe-i,e)| > lACln^^), (10.9) 
16*0*) < expO(l)[^Fft^(^)] = expO(lnln^). 
Using l|10.9f) and the elementary arguments recalled in the second part of Section we get 

hi t d " =omMn jk- (1(U0) 



On the other hand, 



' "" 1} dfl = -r—r{<i)v{k-t){ll>k,s) - </>v(k-l)(Vk-l,e))- (10.11) 



Re 2^Xv#^= ( 10 - 12 ) 



2-ni J lu a v (k-\) 

Combining the real parts of ifTUTjl . ifHrgjl . I|10.10f) (|10.11|) . we get 

1 f G' 
2iri J7 G 



Re(0 y ( M )( j UM, e ) +Sfe=o 1 ( < ? !, Kfe)(/ ife + 1 . s ) ~ 0KM-i)(Mfc+i,s)) ~ ^(o)(mo,s)) 



with the remainder improving to 0(1) if we do not encounter B e . Now, /iM.e = fi>M,e, A*o,s = Mo,s> 
and 

In In 1 

Hk >e = fj.k+i,s + O (h- — (7-r-)), 
m {(j,) h 

with the last remainder improving to 0(h/ In 7^x^) ; if we avoid B e , and (|10.3fl follows. □ 

We end this section by some rough estimates on the location of the skeleton and the corre- 
sponding distribution of eigenvalues for the reduced operators constructed in Sections El El Our 
starting point is the reduced symbol in H2.6(l and the corresponding 1-dimensional symbol 

Q (r,x,£,e, —;h) ={q)(T,x,0+O(e) + — P2 {t,x,0 + hfr + h 2 p 2 + ... . (10.13) 



Here we shall take r real (and eventually of the form h(k — ^p) — f 2 -)). If z is the original spectral 
parameter, we introduce the new spectral parameter w, by 

z = g(r)+iew, (10.14) 
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and we will work under the assumption h 2 <C e <C ft, 1 / 2 . 

Recall the microlocal normal form for Q near the branch point, given by Proposition 16.21 and 
in particular l|6.13|l : 

U^QU = K £th y e (T,I;h) +£>(e-™), I = -{x o hD x + hD x o x), (10.15) 

where the leading symbol in K^^iu is k €i h?r e (r, i) with i — x£, given in Proposition 16.11 Corre- 
spondingly, we replace w by the new spectral parameter fi, given by 

Ke.,h?/e{T,ll\h) = w. (10.16) 

We next estimate the location of the skeleton in the /i-plane, and start with the case > Ch. 
Assume for simplicity that we are in the case 1: Im/i > — C|Re/i|. We will only be concerned 
with the horizontal part S' of the skeleton. When Re // > 0, it is given by the curves T 3A + = T12, 
r 2i4 + = Ti i3 in Ij8.12jl . where 

X(jj.) = r(^) + |Re M , 

YU) = (Re^)arg(-) - Im/i + /iRe0_ (-). 

i fi 

Clearly X(fi) is uniformly Lipschitz continuous and for fi > 0, we get 

X(n) = hReO-(-). 

M 

According to fity . we have ReO_(|) = 0(er 2 \^/ h ). 

When e = 0, h 2 /e = 0, we know that the leading part of Q in i|10.13|) is real-valued (assuming 
that (q) is real for simplicity), so it follows in this case that when fi is real, then ImSjk — 0{h). 
Since Sj t k depends holomorphically on //, we conclude that in general 

Im £?j,fc(At) = 0{e + h 2 /e), /i£R. (10.17) 

Now combine this with (|8.12l) . the estimate X([i) — 0(/ie _27r '' i '/' 1 ) and Proposition lB.ll to conclude 
that in the region \fj,\ > Ch, Re /j, > 0, the horizontal part S' of the spectrum is given by the union 
of two curves of the form 

lmfi = f(Ren), (10.18) 



with /' satisfying (|B.21|) . and further, 

\f(x)\<c(e+-)u i: vs I -y ' , I . (10.19) 

11 [as| 11 ( e +M) 

In the left half-plane, we recall that S' has a more complicated structure. Assume, to fix the 
ideas, that Re fiA < 0. Then S' is the union of the curves (defined in l|8.13|l ): 

Ti 3 : Im n = ji 3 (Re u) (Im u) In - — - = Im Si 2 + (X — 7rRe fi) + 7rRe /1 

H 

r 3 4- : Im /1 = 73 4 - (Re fj) O- (Im fj) In - — - — Im S3 4 + (X — 7rRe a) — 7rRe u, 

\fi\ 
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in the region Re/i a < Re/i < 0. Here 73,4- (Re /i) < 7i,3(Re /i) and the two curves cross at lla- In 
the region Re /x < Re /ia S' is given by 

r 14 - : Im /1 = 7x,4- (R- e A 1 ) ^ (I m A 1 ) m ] — r = 7: (I m ^1,2 + I 7 ™ £3.4) + (-X" — ^Re /i), 

and this curve also contains /iA- 
When fx < 0, we have 

X - 7rRe M = «teO_(~) = ft,C(e- 27r ^l/ h ), 
A* 

so again / := 71.4- satisfies (|1U.19|I . while 

7 3i4 -(Re/i) := ImS 3)4 + (X - nRefi) < 7 3)4 - (Re /i) < 71,3 (Re/i) 

< IraSx,2 + (X — ?rRe/i) =: 71,3 (Re/i), 

for Re/i a < Re/i < 0, where / = 71,3,73,4- satisfy 1)10.191) . 

In the region |/i| < C7i the horizontal part of the spectrum is a union of curves r^fc given in 
H9.4[l . (|9.5|l . Here the new function X is uniformly Lipschitz and O(h), so the skeleton is here 
contained in a region 

(e+-) 

|Im/i| < 0(1)- t~ — . (10.20) 

ln 7T 

The overall conclusion is that the skeleton is contained in a region 

|Im/t| <0(l)(e+— )max (- ^ — ,= — V" | , (10.21) 

where we recall that (Re/i)^ = (h 2 + (Re /i) 2 )^- 

We end this section by establishing a simplified statement, to be used in Theorem 1.1. We 
shall simply remove a small rectangle around /i = where we have seen that the description of the 
spectrum is more intricate. 

Start by recalling the definition of [Xa, Mb prior to Q8.25JI . For instance [ia is the intersection 
of the curves 

A: — 27rRe/i = Im5x,2 — Im 63,4, 
and 

T 14 -: Im ll = 7 14 -(Re/i), where / = 7x 4- satisfies (I10.19|l . Using that 177151,2,17715*3,4 = 
0(e + h 2 /e + \Imfi\), we get 

f h 2 \ f h 2 \ ( I 1 \ 

Re ijla — O [ e H h |Im lia\ I , Im lla = O I e H I max 



lnlRe/iAll' |ln(e+f )|i ' 



implying, 

7,2 0(e+^) 
e |]n( e+ i.)| 



Re/iA = 0(e+— ), Im/iA= ,, , , h Y NI - (10.22) 
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We have of course the same estimates for \xb- 

Choose C > sufficiently large so that the "black box" 



B = [-a, a] + i[-b, b], with a = C{e + — ), b = C M \ : \' 2 , (10.23) 



I Me 

contains ha, fiB- Then we have 



Proposition 10.2 The number of eigenvalues in B is 0(j- + -)|ln(e + ^-)|. The eigenvalues 
outside B are exponentially close to Ti A - U IY2 U I\,3- More precisely introduce 

= {/i S Ti A - \ B; a\ + a±- =0, Re // < 0}, 
Ei,2 = {ti e Tx.2 \ B; 01 + a 2 = 0, Re /i > 0}, 
#1,3 = {ti 6 r x ,3 \ 6; ai + a 3 = 0, Re fi > 0}. 

TTien £/iere is a bijection b (possibly after a slight modification of B) between the set of eigenvalues 
outside B and E 1A - UE 1:2 UE 1:3 , such that b(ji) - fi = 0{e-^ Re ^/ h h/\ In \fj,\\). 

Proof: We may first notice that we can replace the index 4 _ by 4 without changing the validity of 
the statement of the proposition, since 04 — a 4 - = C(e~ 27r ' Re ^/ h ) 7 Re/i <C —h. In view of (|10.21fl . 
we know that there are no eigenvalues outside B with |Re juj < a and the discussion in Section [H] 
then shows that the eigenvalues outside B have to be exponentially close to T14 U IY2 U I\3 and 
that there is a bijection b as stated. To estimate the number of eigenvalues inside B, we simply 
apply Theorem llO.il with 7 a rectangular contour containing B but contained in 2B and working 
directly with <Z4 instead of a 4 ± . □ 

Consider ^4- of the preceding proposition. In view of Q8.3JI . it is given by the quantization 
condition 

£1,2 + #5,4 + 2/i(ln(-/i) - 1) + ^ + 2«C_(^) = 2n(k + ^)h, keZ (10.24) 

Here we recall from the beginning of Section [S] that the term C_(-) is ~ C\ — + C^-) 2 + as 
^ — > 0. We also know that if a = (e, h 2 /e) denote the small additional parameters in the problem, 
then Sj t k ~ X)o° ^jki^^ ct)h v , for (J, k) = (1, 2), (3, 4), where 5J fe are smooth in a and analytic in 
t. Hence the condition l|10.24|) takes the form 

b 1A -{fi,a;h) = 2n(k+^)h, (10.25) 



where 

61,4- (fx, a; h) ~ V 6? 4 - (/i, a)^, (10.26) 



00 



in the space of bounded holomorphic functions defined in truncated sector: 

Re // < — Ch, \lm fi\ < — (— Re/i), 

with 

b\ 4 _ (/i, a) — 2fj,hx(—fx), b\ 4 _ holomorphic in a full neighborhood of /1 = 0, a = 0, (10.27) 
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and 

% A -(jj,,a) = din 1 -"), v > 2. (10.28) 

Notice that the singularity structure (|10.26f> - HlU.28|) of bi±- is essentially unchanged if we replace 
H by /I = a(/i, a; h)((i + hd(n, a; h)), where a, d are classical symbols of order in h with coefficients 
that are analytic near fi = 0, a = and with a elliptic, Re a > 0, |Im a| <C Re a. 

On the other hand, in the region Re fi < — 1/(7, C 3> 0, we know (and that was done for instance 
in |12j). that the eigenvalues sit on a curve and are given by a Bohr-Sommerfeld condition 

b(n, a; h) = 2n(k + ^)h, k € Z, (10.29) 

where 6 is a classical analytic symbol of order 0: b ~ 6 l/ (/x, a)/! 1 ', and where 

&° (/*,<*) = / (10.30) 

•'Tejrt 

Here 7 cx t(A t , a) denotes a closed loop in the energy surface Q°(/i, a, x, £) — w with w and /i related 
by qiU.ltip . that can be obtained from the real energy curve we get by taking fj, real and putting 
a — 0. Clearly = 61.4- , so our discussion gives detailed description about how one can push the 
standard WKB-construction to the limit \fi\ ^> h in the region Re fj, < 0. 
The same discussion applies to £^1,2, -Ei,3- We get the conditions 

6i ) 2((U, ot\ h) = 2nkh and bi^fi, a; h) = lixkh respectively , (10.31) 

where bj t k, (j,k) = (1,2), (1,3) are defined in the truncated sector Re \i > Ch, \Im /i\ < Ren/C, 
an d b" k h ave the analogous properties to those in l|10.27(l . (|10.28(l . for v > 1, while the first part 
of (|10.27|l should be replaced by the condition that 

bj k (fi, a) — [i In \i is holomorphic near /i = 0, a = 0. (10.32) 

bl 2 is the action along a closed loop inside the appropriate complex energy curve, that can be 
obtained by deformation from the case /i > 0, a — where we take the left real component, close 
to the left loop in the oo-shaped set -Kq.o- F° r 4 w ^ deform from the right real component. 



11 Justification by means of a global Grushin problem 

One dimensional Grushin problems. We may assume here without loss of generality, that (q) is 
real-valued. Then we know that / in (|6.14|l satisfies 

00 

f{w;h)~J2f k {w;h), (11.1) 


where /o is real- valued when e, h 2 /e = 0, r € R. Recall that / and fk depend analytically on r 
and smoothly on e, fr 2 /e. 

Also recall that the spectrum of Q is localized to the region (|6.8|l : |Im kj| = C(/i + e + h 2 /e) (as 
follows also from the more refined estimate and in view of (|6.15Jl : f(w; h) — —fi, it follows 

that fi is localized to a domain of the same type. 



58 



We shall now introduce three different Grushin problems for Q — w in the spirit of JT], |28|. 
Let x S Co°(( ne igh (0, 0), R 2 )) be equal to one near (0,0). We realize x as an /i-pseudodifferential 
operator, that we also denote by %, using a Gaussian resolution of the identity (see JS], 
so that our calculus errors will be exponentially small rather than just 0(h°°). Assume, in order 
to fix the ideas, that the support of V% is a thin annulus around (0,0), containing the points 
a.j, j = 1, ..,4 (see Section EJ) Recall the definition of gj in I|(i.l7|l . Let g* be the corresponding 
functions defined for Q* — w, depending anti-holomorphically on w. Define 

Bi + u={-iy{u\~[Q\x]9* j )wr ( n - 2 ) 

Here Wj is a small neighborhood of ctj and (m|w)via, = (Xj u \ v ) where \j € C§°(Wj) is equal to 1 
near ay and we also let Xj denote the corresponding Gaussian quantization. We may normalize 
the choice of g* so that 

R? +9j = 1. (11.3) 

Expanding the commutator, we see that the definition of R 3 + u (up to an exponentially small error) 
is independent of the choice of x, provided that u is a microfocal null-solution of Q — w. 

Assume for simplicity that x i s real-valued and that the corresponding quantization is selfad- 
joint. Put 

i£u_ = (-!)*«_ i[Q,x]&. (11.4) 



Our first Grushin problem is directly adapted to the derivation of the quantization condition 
lj6~2"gj) in Section El It is the form 

( (Q - w)u + = v (115) 

\ R + u = v+, 

with 

R+ = Rl : L 2 d , C, i?_ = Rl_ : C L%. (11.6) 

Using SectionElwe see as in ^]|2H1 that it is well posed for w in some fixed complex neighborhood 
of with a solution of the form 

( u = Ev + E + v + , 

[u- = E-v + E- + v + . 1 ' 

Here we get 

E_+v+ = ihv + {c 2 ^e 2 ^ e M + i^ +s ^ + c 2A e 2 ^+i s ^ - Cl , 3 e 2 ^ + i s ^ - c M ), (11.8) 

where the parenthesis is the same as in the quantization condition i|(j.29|) . As usual, we read off 

the approximate eigenvalues as the zeros of E |_. 

The drawback with this first Grushin problem is that the solution operator (|11.7|) will grow 
exponentially when fi > 0. This can be seen either directly from the explicit formulae for cij.k and 
the slightly less explicit expression for c^fe, or from the fact that for /i > 0, we have approximately 
a double well problem and with R + we prescribe the solution u in l|11.5fl in one of the wells, 
and hence u will in general be exponentially large in the other well. Of course, we will have to 
accept some exponential growth with a rate 0(\lm w\ +e + h 2 /e) but certainly would like to avoid 
exponential growth with a fixed rate when fi is real. 
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It seems impossible to cover a full neighborhood of w = with a single Grushin problem whose 
solution operator does not exhibit exponential growth in some region, so we shall use 2 Grushin 
problems where one will be nice roughly in the upper half-plane and the other in the lower half 
plane. 

The second Grushin problem is designed to cover the region Im /i > with some margin. It is 
of the form Ijll.Sjl with 

R+u = (R 2 + u,R\u), R-U-=R l _u 1 _+R?_v?_, (11.9) 

so that R+ = Lg, — ► C 2 , i?_ : C 2 — > Lg,. For the corresponding model problem for Pq — /i in 
Section H we get (cf(Q~TTJ), g3S> : 

' Ul ^-[/f U2 V (11.10) 



M 3 / \ "4 



^=i^ 2 ! M> i = r(i "f c-v^v ^wi- (ii.li) 



,"3,2 "3,4/ y27T 

This is basically the approach of |11| and as there, we see by direct calculation or by a more general 
normalization argument that U (/i) is unitary when /i is real. We also see that U is uniformly 
bounded in any disc D(Q, Ch) for 

Ch 

lmn> r . (11.12) 

ln 7I 

For /z outside an angle around — i[0, +oo[ with 3> h, we apply Stirling's formula (|5.2|) and 

get 

Uj , k = exp (*t -*t l n( - iM) ±(lt +i * ) + (11.13) 

with the + sign valid for 1*1,2, "3,4, and the — sign for "1,4, "3,2- As in Section |H1 we get 

Kk\ = exp ( - ~(1 + lnirjlm/* + ^(arg(^) ± I) + ReO_(-)). (11.14) 

We now also assume that < C^ 1 <C 1, so that ln(|/i| -1 ) 3> 1. We shall estimate the exponent 
from above. When Re [i > 0, the worst exponent is the one with +7r/2 in the middle term and we 
approximate 

1 , ,AS . i\ !/„ w s 1 (Re^)(Im/x) 
-Re M (arg(-) + -) = -(Re M )(arg M ) - 5 ^ , 

which is dominated by the first term and hence 



%,fe| <exp(— j^ln— +0(l))lm/x), (11.15) 



when Re/i > 0. When Re < 0, the worst case is the one with — ir/2 in the middle term and we 
approximate 

1 n it -Re/i 1 |Re/x|(Im/x) 
-ReMarg(-) - -) = — (* - arg M ) ~ - , 

leading to l|11.15|) also in this case. We conclude that U is bounded in a domain of the form 

Ch 

{ M ;Im M >-- j—}. (11.16) 

in ; .I 1 
h+\pi\ 
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As a consequence, we get 

Proposition 11.1 The problem (|11.5|) with R± given by (|11.9|) is microlocally well-posed (with 
errors 0( e -i/(C*))) for \n\ small with H in D(Q, j) or away from a small conic neighborhood of 
iR . // we write the solution as in (|11.7|) . then for /i as in ((11.16(1 . we have 

h\\E\\, \\E±l \E_ + \ =0(l)exp^(|lm/i|+e+-). (11.17) 

h e 

Let us also compute E y. Near the branching point, we recall that we have the relation ((6.18(1 

for null- solutions to P — z, equal to Ujgj near ay. To determine E |_, we consider (|11.5() with 

v = 0, so that uj = vf fo r j — 2, 4. We then want to express Uj, j — 1, 3 in terms of u 2 , u A . We 
can do this using . 1U|) . (|1 1 . 1 XI) above and redo some of the work in Section |SJ but it is easier to 
to use the work already done and "solve" ((6.18(1 . We get 



/ Cl - 3 n — Cl >3 c 2,4 
"1 \ [ C 2 ,3 C M C 2 ,3 

1 C 2 ,4 



U3 



C2,3 C 2 ,3 




(11.18) 



Notice that in our Grushin problem (with v — 0) are the discontinuities we obtain at ax, 

«3 when trying to extend a null-solution near (0,0) with prescribed u = v^gj near ctj, j = 2,4, to 
a global null-solution near Kq q. We get 




where 



It follows that 



2TT0! = 2tt9 1 + yS 3 4, 2TT02 = 2TT0 2 + ySi 2 . (11.20) 

h h 
dct£_+ = —(e 2 ^i+02) C23 + C24e 2^1) 2 _ c _ j (1L2 i) 

C2,3 

_ e 27ri(5i+5 2 ) _|_ C 2,4 c 27r»§ 2 _ £l£ e 27riei _ C l,4 
C2,3 C 2j 3 C2,3 

From the middle expression, we see that this determinant is equal to the expression in the quan- 
tization condition (|6.29() times a non-vanishing factor. 

The third problem is designed to cover the region Im \i < 0. It is of the form 111.5(1 with 

R+u = {R\u, R\u), = R 2 _u 2 _ + Rtut. (11.22) 

For the corresponding model problem for Pq — /i in Section 01 we get (cf 1)4.11(1 . ((4. 13(1 ): 

^ =v( Ul ) =U- 1 ( Ul ), (11.23) 



with U as in 1(11.10(1 . unitary for real /z, so that 

V(»)=U(-py. (11.24) 
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Thus from (|11.15f) we see that the matrix elements Vj k satisfies 



Vj,k\ <exp-(ln r - r + 0(l))lrn/i, 
h \fi\ 



(11.25) 



and V is bounded in a domain of the form 



r T Ch i 

{H] Imp, < - 1 }. 

111 h+m 



(11.26) 



Proposition 11.2 The problem (|11.5I) with R± given by (lll.22|l is microlocally well-posed (with 
errors 0( e - 1 ^ Ch ))) for small with /i in D(0, j) or away from a small conic neighborhood of 
iR + . If we write the solution as in ((TTTjl . then for /i as in $11. 26)) . we have the estimate l|11.17[) . 

We now compute the corresponding E |_ so we put v = in l|11.5fl and repeat the arguments 

above. Near the branching point u is a null-solution, = Ujgj near ay, now with vf — Uj, j = 1, 3. 
By "solving" (|6.18|1 . we express U2,«4 in terms of ui,U3 and find 



We then get 



= E. 



C2A 


Cl,3C2,4 \ 




C2.3 




ClA 


Cl,4 




1 


Cl,3 


")< 


Cl,4 


Cl,4 




/ c 2,4 




C2,3 - 


Cl,4 




V 


1 


£M + 


Cl,4 


C2,4 



Cl,3C2,4 

Cl,4 
„— 27T101 



(11.27) 



(11.28) 



and 



dct£L_ 



-2TTi(e 1 +0 2 



ClA 



C2.3 
Cl,4 



c 2,4 27Ttfli 
Cl,4 



Cl,4 



27rie 2 

-27ri(6li+e2; 



27ri.6 



Cl,4) 



(11.29) 



The global Grushin problem. We first explain what the natural range will be for e. Our global 
Grushin problem will be built from a direct sum of problems for the operators P e — z in l|2.4|l . 
These operators can be written 



g( h (k- -r) - tt) + itQ{h{k- -£-),x,hD x ,e, —;h), 
4 2n 4 e 

where Q is the operator appearing earlier in this section and in Section |SJ 
Clearly, we add conditions "R±" only for such k for which 



(11.30) 



Z-g(h(k-!f)-f) 

it 
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is close the spectrum of Q, i.e. for which 

h 2 

\lmw k \ <0(e+— ), (11.31) 

(cf. 1(6. 8() ). assuming for simplicity that (q) is real- valued. Then according to (|11.17fl we can expect 
that our Grushin problem will have an inverse E with norm ||£|| = 0(\)e c ^ +! i\ 

Now we cannot expect to have a complete decomposition into a direct sum of operators (|11.3Ufl , 
but Proposition 13.31 shows that so is possible up to an error C(l)exp (— l/(C(e + h))). Thus in 
order to absorb the error by a standard perturbation argument, we need 

which would follow from 

1 e h , 
—— » - + -, 11.32 

e + h h e 

or equivalently 

e 3 + h 3 < eh, 

Here we already know that h 3 <§C eh, since we assume h 2 /e <C 1, so the new constraint is e 2 <C h. 
From now on we work in the range 

/i 2 <e</i 1/2 . (11.33) 

Since g is real-valued with g' ^ 0, it follows from l|6.8(l (or 1(6.9(1 in the general case, when 
(q) is not assumed to be real) that the operators (|1 1.30(1 have disjoint spectra. When 1(11.31(1 is 
never satisfied for any k, it is straight forward to see that z is not in the spectrum of our original 

operator P c . Assume now that (|11.31fl holds for (at most) one k = k. Let be the corresponding 
operators R± defined earlier in this section. Using the notation of Proposition 12. II we define 

R+u = Rfde^U^e-^uler^j^)^^)) (11.34) 
R-U- = e^ G Ue'^ A {er hSL _s^®R { ^ ) U-), 



2nh 



for u £ L 2 (M), U- G C. Here e k (t) = e lkt . 

Repeating the arguments from ^1 ^] we see that we get a well-posed problem with 

E_ + (z)=E^l(z) + 0(evr£v). (11.35) 

(k) 

Here E_! is the E \. of the approximate 1-dimensional Grushin problem treated in the appro- 
priate one of the Propositions 111.11 111.21 depending on the sign of the corresponding parameter 
Im fi. 

In the expressions <|11.21|) . 111.29(1 . we have the term e ± 27ri ( 6 'i+ 9 2) i wne re 9j are given by ((11.20(1 . 
We see that Imfy = 0(f + |), so 

e ±!hri(? 1 +? a ) =e O(t+i). (11.36) 

We conclude from this and 1(11. 32|) that the remainder in 1(11.35(1 is 0(eT cte+M ) times the domi- 
nating term in the expression for £'i fc | in 1(11. 21|) , or ((11.29|) respectively. This implies that if we 

pass to the [i-variahle (for k = k) and define the skeleton as in Sections^ and the corresponding 

body B as in the beginning of Section \l(A then the zeros of det E |_ are confined to B and Theorem 

still applies to give the number of eigenvalues (in the /i-plane) inside an admissible curve. 
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12 Improved parameter range for barrier top resonances in 
the resonant case 



We start the discussion in this section with the following general observation. Let P e be a smooth 
family of operators, satisfying all the assumptions of the introduction, and in particular l|1.13fl . 
In Proposition 2.1 we have seen that the operator P e= o can be reduced by successive averaging 
procedures to 

P e =o = g{hD t ) + hpi(hD t , x, hD x ) + h 2 p 2 {hD u x, hD x ) + ...., t e S\ r, x, £ « 0. (12.1) 

Proposition 12.1 Assume that the subprincipal symbol of P e =o vanishes and that the spectrum 
of P t= o clusters into bands of size < 0(l)h N ° , for some integer N > 2. Then Pj{r, x, £) = Pj{r) 
are independent of (x 1 £) for 1 < j < N — 1 in (|12.1I) . 

Proof: Since the subprincipal symbol vanishes, we already know that p\ = 0. Suppose that the 
conclusion of the proposition does not hold and let Ni E {2,3, N — 1} be the smallest N 
with Pn{ t 0i x, £) non-constant for some r ~ 0. Take a family of Gaussian quasimodes e a (x), 
a = (a x , a^) € neigh (0, R 2 ) with 

||e Q || = 1, p'% 1 (ro,x,hD x )e a = p Nl (tq, a)(a)e a + 0(h 1/2 ) in L 2 . 
See, e.g. |S] for the standard construction of such quasimodes. Then put 

f a , h = {2^ 2 eM h ( k -^)-^>e a {xl 
with k — k(h) such that h(k(h) — feo/4) — Sq/2tt — > To, so that 

KU, k =9 (h (*" " i) +/lAri ^ (* (*" T) ~ I'") + in ( 12 - 2 ) 

Hence, since we are dealing with selfadjoint operators, 

dist [g (h + (fc (* ~ ~ «) , < 0(^ 1+ *), (12.3) 

and varying a, so that the values pjq t (h(k — ^Q.) — j 2 -, a) fill up a whole interval, we get a contra- 
diction to the clustering assumption. □ 



Remark. Proposition 12.1 remains to hold in the case described in Section 4 of where the 
operator P e= o is conjugated into a normal form in a neighborhood of a Lagrangian torus, rather 
than near a closed iZp-trajectory. 

From now on we shall assume that P e =o satisfies the assumptions of Proposition 12.1. Let us 
now switch on e. An application of Proposition 2.1 together with Proposition 12.1 then shows that 
microlocally, near a closed Hp-trajectory, P e can be reduced to the form 

P t = g(hD t ) + e (i(q)(hD t ,x, hD x ) + 0(e) + O ( —) + hpx + ..) . (12.4) 
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It follows therefore that in the results of ^21 (Theorems 6.4 and 6.7 there), we can replace the 
exponent 2 by the exponent Nq in the parameter range for e. Thus for the study of the spectrum 
of P e in a region where \Ke z\ < 1/(9(1) and \lmz/e — Fo\ < 1/0(1), when Fo is a non-critical 
value or a non-degenerate maximum or minimum of Re (q) along p _1 (0), it suffices to assume that 

h N ° < e < h 5 , (12.5) 

for some S > 0. In the case when Fo is a saddle point value of Ke(q), from Theorem 1.1 we get 
the condition 

h N ° < e < h 1/2 . (12.6) 

Indeed, in the latter case we still have Proposition 13 . 31 and the decoupling condition analogous to 
<|11.32[l becomes e 3 + h No+1 < eh, which is fulfilled by (|l2~B|l . 

We shall now apply these observations to improve the result of Proposition 7.1 of ^2]j giving 
a description of the individual barrier top resonances of the semiclassical Schrodinger operator in 
the resonant case. Before doing so, and also for the future use in Section 13, we shall first briefly 
recall the general setup in Section 7 of ^21> as wen as m Section 5 of [T3| . 

As in ^1 EH , let us consider 

P = -h 2 A + V(x), P{x, = £ 2 + V(x), (x, e T*R 2 , (12.7) 

where V satisfies the general assumptions of Section 7 in ^2] > allowing us to define the resonances 
of P in the lower half-plane inside some fixed neighborhood of Eq > 0, where V(0) — E$, V'(Q) = 0, 
V"{0) < 0. As in H21 ESI, we assume that {(0,0)} is the only trapped Hp-trajectory in P- 1 (E ), 
After a linear symplectic change of coordinates, we may write 

2 A 

P(x,O-E = Y / Y^~ x ^ +P3 ^ +p ^ + -^ foO-^O, (12.8) 

where Xj > and Pj(x) is a homogeneous polynomial of degree j > 3. Recall further from |12| 
that the study of resonances of P near Eq can be reduced to an eigenvalue problem for P after 
applying some variant of the method of complex scaling, and that near x = this simply amounts 
to working in the new real coordinates (a?, £), given by x = e l7r / 4 x, £ = e~ J7r / 4 £. 

Performing the scaling and dropping the tildes from the notation, we see that the problem 
reduces to studying the eigenvalues close to of the operator i{P — E ), now elliptic outside a 
small neighborhood of (0, 0), with symbol 

P(x, = p(x, + ie^pz{x) + ie" Pi {x) + ie^p b {x) + (12.9) 

where 

P(^) = Ey(^^ (12-10) 

Here we continue to write P to denote the scaled operator. 
We assume that Xj > in (|12.9|) are rationally dependent, 

3fc° = (fcj, fc°) e Z 2 \ {0}, fcjAi + k° 2 X 2 = 0, (12.11) 

which implies that the iJ p -flow is periodic. 
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As in |13| . we are interested in eigenvalues E of P with \E\ ^ e 2 : Q < e ^ 1. After a 
rescaling x — ex, and dropping the tildes over the new variables, we get an operator P e — j?P that 
we view as an /i-pseudodifferential operator with the symbol 

P e (x, = \p{e{x, 0) = p(x, + iee 3 ^p 3 (x) - ie 2 Pi (x) + 0(e 3 ). 

Here h = h/e 2 . Now the spectrum of P e =o is that of the harmonic oscillator, and hence it clusters 
into sets of diameter and separation of order h. An application of Proposition 12.1 shows that 
all the pj in l|12.1(l are constant. Moreover, since in this case all the eigenvalues depend linearly 
on h, we see from the proof that the pj have to vanish. It follows from l|12.5|) that in the zone 
corresponding to non-critical values Fq or non-degenerate maxima or minima, the range of energies 
that we get is 

h ^ a /(l+2N ) <<; ^ _ £ Q | < /l 25/(l+25) ) (12 12) 

for all N — 2, 3, . . . and all 6 > 0. When F corresponds to a branching level, we get from (|12.6|) 



h 2N /(l+2N ) KlE-Eol^h 1 



/2 



(12.13) 



We summarize the discussion above in the following proposition, which is an improvement of 
Proposition 7.3 in |12j . Clearly, in a similar fashion, we also obtain an improvement of Theorem 

5.1 in nni. 

Proposition 12.2 Assume that the principal symbol P(x,£) in (|12.7|) has an asymptotic expansion 
l|12.8f) . and assume that (|12.11|) holds. Assume furthermore that the function (pa), defined as 
the average of p 3 along the Hamilton flow of p in <|12.10fl does not vanish identically. Then the 
resonances of the operator P in the domain 



{z e C; h 2No/{1+2No) < \z-Eo\ < ^}\IJ{ Z e C; Rez " E Q -A\lmz 
where rj > 0, S > 0, and No = 2,3,... are arbitrary but fixed, are given by 



3/2 



< r\ \Im 2| 3 ^ 2 | , 

(12.14) 



Eo-il h(ki - a/4) + e 3 J2 W e 



with 



2'., 



3=0 



4u 



r =ie 3 ™/ 4 (p 3 )(0+O[e + 



e l+2AT 



k 



h N ° 

d+2N a 



r 3 =0(l), j>l. 



(12.15) 



We have k = (fci,/c2) € Z , S — (Si,^) with Si = 2n, and a = (01,0:2) G Z 2 is fixed, and we 
choose e > with \E — Eq\ ~ e 2 . The union in (|12.14|) is taken over the set of critical values of 
(ps), restricted to p _1 (l), with A varying over this set. 



Remark. If (p^) restricted to p _1 (l) has precisely one non-degenerate saddle point with the critical 
value A, then the results of the present paper apply and give a description of the individual 
resonances in a half-cubic neighborhood of the curve Rez = Eq + A\Im z\ 3 ^ 2 . In the following 
section, we shall consider explicit examples of homogeneous polynomials for which the assumptions 
of Theorem 1.1 are satisfied. 
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13 Examples in the barrier top case 

This section is devoted to a study of examples of potentials of the Schrodinger operator (|12.7|) to 
which Theorem 1.1 is applicable. 

Let us recall from l|12.9fl that we are interested in eigenvalues close to of the operator P, 
elliptic outside a small neighborhood of (0, 0) with symbol 

P(x, £) = p(x, + ie^p 3 (x) + ie™ Pi {x) + ..., (13.1) 

where the harmonic oscillator p(x, £) has been defined in (|12.10(l . As before, we make the resonant 
assumption (|12.11|) . 

Consider first a general perturbation of p of the form of a linear combination of terms x a t;° 
with \a\ + \j3\ = to, for some m £ {3,4,5, ...}. Recall from ^21 how to compute the corresponding 
trajectory average Basically we use action-angle coordinates, but to start with, we can do 

things a little easier by introducing 

Zj = Xj + i£j e C, (13.2) 
and notice that along a _ff p -trajectory we get in the Z\, z 2 coordinates: 

Zj(t) = e- iX ^Zj(0). (13.3) 
Then write Xj (t) = Re Zj (t) , £j (t) = \m, Zj (t) , so that 

2 

x(t) a £{t)P = Y[((Re Zj (t)) a Hl m z 3 (t))^) (13.4) 

i=i 

2 



l ww n((z,(0) e -^* + z7(0)e^ t )^(z,(0) e -^ t - 2-(0)e^*)ft) 



2l«l + l/3UI/3| 

Then expand the product by means of the binomial theorem. The time average is equal to the 
time-independent term and since this average is constant along each trajectory we shall replace 
the symbols Zj(0) simply by Zj. 

In this section we consider the case when A x = A 2 = 1, m = 4, = 0. (In we noticed 
that in this case the average will vanish when m — 3 and in |13) we made a more refined study 
of that case taking into account one more term in the perturbative expansion). This means that 
we take ps(x) = 0, and for simplicity we also assume that p m — for all odd m in (|13.1(l . so that 
we can concentrate on the perturbation — ip^ in H13.1JI . Performing a rescaling as described in the 
previous section, with e replaced by e 1 / 2 , (i.e. setting x = e 1 / 2 ^ rather than x — ex), we get 

p(x, £) - iepi{x) - e 2 p 6 {x) + ie 3 p s {x) + e 4 p w (x) + (13.5) 

and here, as before, we choose e of the same order of magnitude as the modulus of the eigenvalues 
for the operator P(x, hD), that we want to study. 

Now we continue the calculations of trajectory averages using l|13.4|l . We have 

(xf) = ^((z 1 +z 1 ) i ) = ^(z i 1 +iz 3 1 z 1 +6z 2 1 z 2 1 +4z 1 zl + zt) (13.6) 
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In the same way, we get 



(4) = ?M 4 . (13.7) 



Next look at the averages of mixed terms: 

(x\x 2 ) = ^((z 1 +z 1 f(z 2 + z 2 )) = ^((zf + 3zfz 1 +3z 1 zl+^)(z 2 + z 2 )) (13.8) 

= -^-(\z 1 \ 2 z 1 z 2 + |zi| 2 ziz 2 ) = ^|zi| 2 Re(ziz 2 ). 
lb o 

( Xl xl) = ||z 2 | 2 Re(z 2 zi), (13.9) 

(x\xl) = ^{{zl + 2z 1 z 1 +il){z 2 2 +2z 2 z 2 + zl)) (13.10) 
- ^(^+4| Zl | 2 |z 2 | 2 + z 2 z 2 2 ) 
= lRe(^) + i N 2 N 2 . 

Notice that our averages are invariant under the anti-symplectic involution 

j:{x,0~{x,-0. (13.11) 

This is necessarily the case since we stay in the framework of ordinary Schrodinger operators 
(without magnetic fields) whose symbols have this invariance. 

Now write our results in the action angle variables (pj ,6j), given by 

zj = y/2p]e- ie * , (13.12) 

so that 

(4) = \pl (4) = \pl (13-13) 

{4x2) = \pTpT cos(0x - 6 2 ), ( Xl x 3 2 ) = \pTpT cos(0 2 - 1 ), (13.14) 

(xjxl) = PlP2 + ^ PlP2 cos2{e 1 - 9 2 ). (13.15) 

It follows from the Hamilton equations that pj and 9 := 6\ — 9 2 are constant along every H p - 
trajectory. The involution j can also be described as (zi, z 2 ) (zi,z 2 ), and hence in action-angle 
variables as 

(pi,p 2 ,9 1 ,9 2 ) » {pi, P2 ,-9 1 ,-9 2 ) (13.16) 
We shall study our averages as functions on the abstract symplectic manifold 

£=p _1 (l)/expRHp. (13.17) 
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Using the (zi, z 2 )-coordinates, we have 

p-\l): i(|z 1 | 2 + |z 2 | 2 ) = l, (13.18) 
and the equivalence relation induced by the -Hp-flow is: (z\, z 2 ) <~ (wi,w 2 ) if and only if 

(101,102) = {e it z 1 ,e it z 2 ) 

for some t E R. Thus we see that £ can be identified with the complex projective space P(C 2 ). 
It is well known that this space is diffeomorphic to S 2 . Indeed, P(C 2 ) can be identified with the 
1-point compactification CU {00} via the map (zi, z 2 ) 1— > z\jz 2 and the one point compactification 
can be identified with the Ricmann sphere. 

£ can be parametrized by (pi, p 2 , 9) with pj > 0, p\ + pi = 1, 9 € R/27rZ, with the convention 
that all the (1, 0, 9) denote the same point and similarly for (0, 1, 9). The involution j induces the 
anti-symplectic involution 

j: £9(pi,p 2 ,0)~(pi,p2,-0). (13.19) 

Notice that the set of fixed points of j is given by all points with 9 = or 9 = tt. These points form 
a (great) circle on X and can also be described as the set of trajectories in p _1 (l) whose x-space 
projections hit the boundary of the potential well {x <G R 2 ; \x\ = 1}. 
We consider perturbations of the form 

2 2 
q(x) = -a{x\ + x\) + bx\x\ + -c{x\x 2 + Xix\). (13.20) 

Then on S we get with p = pi, so that p 2 = 1 — p: 

(q) = a{p\ + p 2 ) + b PlP2 + ^pip 2 cos(20) + c{p\ pi +pfpl) cos 9 (13.21) 

= a(p 2 + (1 - pf) + bp(l - + I cos(20)) + c(p5(l - p)l + pl(l - p)5) cos0 

= a+ (6 - 2a)p(l - p) + ^p(l - p)cos(26>) + cp^(l - p)i cosfl 

= a+ - 2a)p(l - p) + fep(l - p)cos 2 + cp^(l - p)3 C os6», 

where we used that p^(l — p)i +pi(l — p)i = p^(l — p)i . 

We are interested in the critical points of this function on S, and the values p — 0, 1 will have 
to be treated separately. In particular we are interested in the number of saddle points. If we have 
only one saddle point we will be able to apply the results of this paper. This is still the case if 
there are two saddle points provided that the corresponding critical values are different. We will 
also encounter the case of two saddle points S\ , S 2 away from the equator and then necessarily 
with j(Si) = S 2 . In that case the critical vales will be equal and the results of this paper will not 
apply directly. We plan to return to that case in a future paper, where the role of symmetries will 
be studied. 

Put 

d=--2a, g = p^(l - p)5, y = C os9. 
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Then, 



(q) = a + dg 2 + bg 2 y 2 + cgy. 



(13.22) 



Notice that y = cos 6* is critical precisely when 8 — 0,tt and that y E [—1,1]. When y ^ ±1, 
we may treat g as an independent variable. The same observation is valid for g(p) 6]0, b]. It is 
non-critical in [0, i[ (i.e. for p ^ ~). (As already mentioned, the value g = 0, corresponding to 
p = 0, 1, will require a different treatment.) 

In order to avoid various degenerations, we shall assume 

d^O, (13.23) 

When c ^ 0, we have b^O, b + d^Q. (13.24) 



1) Critical points with 

0^0,77,^0,^,1. (13.25) 
Here both y and g can be treated as independent variables and the critical points are determined 

by 

J 2bg 2 y + cg = 0, 

\ 2dg + 26. g2 / 2 + cy = 0. 

This can also be written 

f 5 (26. 9 y + c) = 0, 

\ 2d 5 + y(2bgy + c) = 0. 

Under the assumption Ijl3.25|l we have g ^ 0, so we get 

f 26. 9 y + c = 0, 
\ 2dg = 0. 

This is in contradiction with the assumption p3.23|l . so we conclude that there are no critical 
points away from the union of the "vertical circle" given by 8 6 {0,7r} and the horizontal circle: 

2) Critical points on the horizontal circle away from the vertical one: 

0^O,7T, p=i. (13.26) 

Then g = 1/2 and this is a critical value, so we only have to look for critical points with respect 
to y, leading to 

V = -\ (13-27) 
Recall that \y\ < 1 under the assumption i|13.26|) . so we reach the conclusion that if 

l\ < 1, (13.28) 
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then there are two distinct critical points in the region (|13.2(i[) . given by 

P=\, cos0 = -^, (13.29) 

and otherwise there are no such points. (In the remaining degenerate case b = c = the whole 
horizontal circle is critical.) 

We also study the nature of the critical points, by computing the hessian of (q) with respect to 
p, y. Using that g'{\) = 0, g"(h) = — 2, we get at both points 

d 2 y (o) = \, d y d p (q) = 0, d 2 p (q) = -2d. (13.30) 

So both critical points arc of signature (b, —d) where the first component corresponds to the 
horizontal (8) direction. (We use the convention that a signature described by (a, f3) is given by 
(sign (a), sign (/?)).) 

3) Critical points on the vertical circle away from the horizontal one and from the poles p = 0, 1: 

#e{0,7r}, P ?{0, i l}. (13.31) 

For 9 — 0, we have y = 1 and we look for critical points of g i— ► (d + b)g 2 + eg, leading to 

"—why. ^ 

Hence we get two critical points in this region if 

- 1 < th < °- (13 - 33) 

and otherwise no point on this half of the vertical circle. (In the degenerate case b + d — 0, c = 
the whole vertical circle is critical.) 

For 6 = 7T, we have y = — 1 and we look for critical points of g i— > (d + &)g 2 — eg, leading to 



8 = W+d>> (13 - 34) 



so we get 2 critical points in this case if 



< — — y < 1, (13.35) 
b + d y ' 

and otherwise no critical points on this half of the vertical circle. We will see shortly that we have 
critical points at the poles when c = 0. 

In both subcases, we get by a straight forward calculation: 

{q)" gg = 2(6 + d), {q% e = 0, {q)% = 2(fe ^ d)2 4 

so the signature is 

(d + b,d) (13.36) 
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where the first component refers to the direction of the vertical circle through the critical point. 

4) The two points of intersection of the two circles: Here p = | and 9 <E {0,7r}. Here both g and 
y are critical, so our intersection points are both critical. By straight forward calculation, we get 
for C f : 9 = 0, p = \ 

(Q)ee = -l(b + c), (q)'j p = 0, \{q)" pp = -c-b- d, (13.37) 



and for C(, : 9 = n, p = \ 



(q)ee = l(-b + c), = 0, \{q)" pp = c - b - d. (13.38) 



In particular, the signature is 



(— c — b — d, —b — c), when 9 = 0, 
(c — b — d, c — 6), when 9 = n. 



(13.39) 



Also notice that 

{q )( Cf )=a+^ + ^ {q){ C b )=a+^- C - (13.40) 

5) It remains to study the "poles", given by p = 0, 1. Here the p, 9 coordinates degenerate and we 
return to the z-coordinates. Using that 

Re(zfzl) = 2(Rez 1 z 2 ) 2 - M 2 |z 2 | 2 , 

we get 

<<?) = l(\zi\ 2 + N 2 ) 2 + ( b - - l)\zi\ 2 \z 2 \ 2 + b -(Re Zl ! 2 f + ICI^il 2 + N 2 )Re ( Zl z 2 ). 

Make the change of variables Q = Zj/V2 and restrict to the energy surface p _1 (l), which now 
becomes |Ci| 2 + IC2I 2 = 1- Then we get 

(q) = a + d|Ci| 2 |C 2 | 2 + b(Re (C1C2)) 2 + cRe (C1C2), (13-41) 

again with d = | — 2a. 

Recall that we work on the projective space, described as the 3-sphere |(i| 2 + IC2I 2 = 1 modulo 
the action of the rotations t 1— » (e lt (1 , e lt (2) ■ Consider the case p — 0. Correspondingly, we can 
choose the point (CijC^) = (0) !)■ The Hp-integral curve through that point is t 1— » (0, e~ lt ) and 
locally, we can identify S with the transversal hypersurface H in the 3-sphere which is given by 
Im C2 = 0. Thus (2 = 1 — w with w G neigh (0, R), and we get 

w = i-(i-|Ci| 2 )^ = i|Ci| 2 + o(ICil 4 )- 

We can use the real and imaginary parts of £i as local coordinates on H. Then on H , we get the 
Taylor expansion: 

(q)=a + dldl 2 + b(Re Ci) 2 + cRe (i + 0(\(i\ 3 ). 
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We conclude that the "pole" p — is a critical point iff c = and when this point is critical, 
the signature is (d + b,d), where the first component corresponds to the direction of the (vertical) 
circle through the pole. By symmetry in the indices f , 2, we have the identical conclusion for the 
opposite pole, given by p = 1. Notice finally that this case together with the case 3 give a complete 
description of the critical points on the vertical circle away from the crossings with the horizontal 
one. 

We observe that the critical points away from the intersection of the two circles are non- 
degenerate and keep constant signatures under small perturbations of the parameters (except in 
the degenerate cases c = b = and c = b + d = 0) . These critical points can only be killed or born 
by passing through one of the two crossing points. This happens in the following 4 cases: 

1) | = —1: The critical points on the horizontal circle coalesce into the crossing point = 0, p = ^. 
When c/b goes from — f + e to —1 — e, the two critical points disappear and the signature at 
9 = 0, p = \ goes from (— d, —be) to (— d, be) 

2) | = 1. The two critical points on the horizontal great circle coalesce into the crossing point 
9 = 7r, p = i. When | goes from 1 — e to 1 + e, the signature of that crossing point goes from 
(— d, —be) to (— d, be). 

3) = — 1: The two critical points on the vertical circle coalesce into the crossing point 9 = 0, 
p =\- When goes from — 1 + eto — 1 — e, the signature of that crossing point goes from (ce, d) 
to {—ce, d). 

4) = 1: The two critical points on the vertical circle coalesce into the crossing point 9 = it, p = 
\. When goes from 1 — e to 1 + e, the signature of that crossing point goes from (—ce, d) to 
(ce,d). 

In the following, we may assume in order to fix the ideas, that d > 0. In the b, c-plane, we 
define the following open sets, separated from each other by the 4 lines c = ±6, c = ±(b + d), where 
all the critical points will be non-dcgcncratc: 

A: b > 0, -b < c < b. 

B + : max(6, —b) < c < b + d. 

: -(b + d) < c < min(6, -b). 

C+ : c > max(i) + d, — b). 

C- : c < min(fe, —b — d). 

D : b < 0, max(fe, -b - d) < c < min(-6, b + d). 
E + : max(6 + d,-b- d) < c < -b. 
E- : b < c < min(-fr - d,b + d). 
F : b < -d, b + d < c < -b - d. 

Then the earlier discussion gives the location and the signature of the critical points in each of 
the cases. Let C/ denote the "forward" crossing point of the two circles, given by p = \, 9 = 0. 
Similarly let Cb denote the "backward" crossing point, given by p = ^, 6 = ir. 

A: 

Signature at Cf. (— , — ) 
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c = b + d 




Figure 6: When the parameters are in the regions B + , £>_, E + , and E_, we have precisely one 
saddle point, and therefore the results of the present paper apply. In the cases C + and C_ there 
are no saddle points at all, while in the case D there two saddle points. The corresponding critical 
values of (q) are separated provided that we assume that c ^ 0. 
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Signature at Cf,: (— , — ) 
Away from the crossings: 

On the horizontal circle: Two critical points with signature (+, — ) 
On the vertical circle: Two critical points with signature (+, +) 

B+: 

Signature at Cf. (— ,— ) 

Signature at Cf,: (— ,+) 

Away from the crossings: 

On the horizontal circle: No critical points 

On the vertical circle: Two critical points with signature (+, +) 

Here (q)(Cb) is smaller than (q)(Cf) but larger than the two other critical values. 

B_: 

Signature at Cf. (— , +) 

Signature at C(,: (— , — ) 

Away from the crossings: 

On the horizontal circle: No critical points 

On the vertical circle: Two critical points with signature (+, +). 

Here (q)(Cf) is smaller than (q)(Cb) but larger than the two other critical values. 

C+: 

Signature at Cf: (— , — ) 

Signature at C,: (+, +) 

Away from the crossings: 

On the horizontal circle: No critical points 

On the vertical circle: No critical points 

CI: 

Signature at Cf. (+, +) 

Signature at C&: (— , — ) 

Away from the crossings: 

On the horizontal circle: No critical points 

On the vertical circle: No critical points 

D: 

Signature at Cf. (— , +) 
Signature at Cb'. (— ,+) 
Away from the crossings: 

On the horizontal circle: Two critical points with signature (— , — ) 
On the vertical circle: Two critical points with signature (+, +) 

Here (q)(Cf), (q)(Cb) are larger than the values at the critical points on the vertical circle and 
smaller than the values at the critical points on the horizontal circle. From 1)13.40(1 we also know 
that (q)(C f )- (q)(C b ) = c. 

E+: 

Signature at Cf. (— , +) 
Signature at Cf,: (+, +) 
Away from the crossings: 

On the horizontal circle: Two critical points with signature (—,—). 
On the vertical circle: No critical points. 
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In this case (g)(C/) is larger than (q)(Cb) but smaller than the two other critical values. 

Signature at Cf. (+, +) 
Signature at Cb'. (— ,+) 
Away from the crossings: 

On the horizontal circle: Two critical points with signature (—,—). 
On the vertical circle: No critical points. 

In this case (q)(Cb) is larger than (q)(Cf) but smaller than the two other critical values. 
F: 

Signature at Cf. (+,+) 
Signature at Cb- (+,+) 
Away from the crossings: 

On the horizontal circle: Two critical points with signature (—,—). 
On the vertical circle: Two critical points with signature (— , +). 

In the cases B + , £>_ , E+ , E- we have precisely one saddle point (necessarily) situated on the 
vertical circle which is the fixed point set of j. In these cases, the results of this paper apply. The 
results also apply in the case D, provided that we assume that c ^ in order to separate the two 
saddle point values. In these cases it is easy to understand the structure and the shape of the level 
sets (q) — C. In particular we see that when we let C be a saddle point value we get a connected 
"oo" shaped set (and no "circular" components on which (q) is non-critical everywhere). 

In the case F, we have two saddle points situated on the vertical circle symmetrically with 
respect to the horizontal circle. Since we have chosen to use perturbations which are symmetric 
under permutation of xi,x 2 , the function (q) is invariant under the map p ► 1 — p, so the critical 
values are necessarily equal. Here we can break the symmetry by adding a small multiple of for 
instance x\ so that we still have precisely two saddle points, but with different critical values. Then 
the results of our paper apply. 

In the case A, we have two saddle points on the horizontal circle. They are of course exchanged 
by application of j and this symmetry remains under perturbations within the class of Schrodinger 
operators without magnetic field. We hope to analyze this case in a future work. 

We shall next compute the e 2 -contribution to the averaging of the principal symbol p(x, £) — 
iepi(x) — e 2 p 6 (x) + 0(e 3 ) appearing in in l|13.5[) . by applying the calculations of the end of Section 
El with 

q = -p±{x), r = -p 6 (x), T = 2tt. (13.42) 
Recall from there, that we have the averaged symbol 

Pe| A ^ P + ie(q)+e 2 ({r) - ~C(q,q)j +0(e 3 ), (13.43) 

where C(qi,q 2 ) and Cor(qi,q 2 ) were defined in (|8.ti9|) . I|8.67|l . 
A simple calculation gives 

{zV 3 }, = 0, {z a ,z?} =2i(^± ^\ ^ 

More generally, 



N 2 M 2 



r a -5 B-0 X _ 9 . / g(ai,ai;/3i,/3i) afa, « 2 ; flz, fo) \ ~+/j~s+0 

1,1 |z!i 2 + |z 2 p r ' 
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where a denotes the symplectic form, viewed as an alternate bilinear form on T*R 2 x T*R 2 . Hence 
Cor(^,^; S ) =2i( l^h^hM + ^p^p^m-M) 

\ N N 2 J 

when \a\ - \a\ = |/3| - |/?|, and Cor (z Q z 5 , z%^; s) = otherwise. If a = (ai, a 2 ), 6 = (&i, 62) € N 2 
with I a. I = |6| =4 we get by multinomial expansion 

Cor (x a ,x h ) = -g Cor ((2 + 2)°, (z + z) b - s) (13.44) 



1+5=0 

/3+/3=6 



2/ 

28 



& v a / v^/ M 2 N 2 ; 



c + 5 = 
/3+j8=b 
-|c| = |/3|-|/3 



When calculating this kind of expressions, it is useful to observe that the relations \a\ + \a\ = 
\f3\ + \0\ = 4, |5| - |a| = \0\ - \(3\ imply: \0\ = \a\, \0\ = \a\. 
We have the Fourier series expansion 

k£Z\{0} 

Combining this with (|13.44l) and the Parseval identity, we get 
1 f 2lT 

C{x a ,x b ) = — / (s-7r)Cor(a; a ,a; b ;s) (13.45) 
27r Jo 



a\ ( b 



2_ x - \aJ\Pj I a(a 1 ,a 1 ;P 1 ,(3 1 ) <r(a 2 , a 2 ; A, AO ) a+p-a+p 

^ I5I-H U1I 2 + U2I 2 



z ' z 



a + a = a 
f> + fi = b 
|5|-|"I = I0|-|/§I 

|a|-|a|^o 

Using this formula we get after a few days of simple but tedious calculations: 

C(xi + x 4 2 , x\ + x\) = -^(Nil 6 + N 6 ), (13.46) 

C{x\ + xi,xlxj) = -^(3\z\ 2 (z 2 z 2 + z 2 z 2 ) + 16\ Zl \ 2 \z 2 \ 2 ), (13.47) 
C{x\+xi,x\x 2 +xxxl) = ^(2(^3 + z3 z 3)_(51(|zi| 4 + |z 2 | 4 )+36|zi| 2 |z 2 | 2 )(ziz 2 + ziz 2 )), (13.48) 

C(x\x\,x\x\) = -l|z| 2 (9|z!| 2 |z 2 | 2 + 8(z 2 z 2 + z 2 z 2 )), (13.49) 
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C(x\x 2 ,x\x 2 + x x xl) = -t^((17(|zi| 4 + |z 2 | 4 ) + 90|zi| 2 |z 2 | 2 )(z 1 z 2 + z x z 2 ) + 12(2? *| + zf4)), 

(13.50) 

C{x\x 2 + Xl xl x\x 2 + Xl x\) = -l(17(|z!| 6 + |z 2 | 6 ) + 153\ Zl \ 2 \z 2 \ 2 \z\ 2 + 51\z\ 2 (z 2 z 2 + z\z 2 )). 

(13.51) 

Now recall that q is given by (|13.2Ufl . so that by l|13.43fl . we have 

Pe\ AeG p + ie{{q)+ief(r,a,b,c) + 0(e 2 ))=:p + ieq e , (13.52) 

1 A 



f(r,a,b,c) = -(r) + ^qa 2 C(x 4 1 +x 4 2 ,x 4 1 +x 4 2 ) + b 2 C(x 2 xlx 2 x 2 ) + 

^-c 2 C(xlx 2 + xixl, x\x 2 + xixl) + ^-C(x 4 + x\, x\x 2 2 ) + 
9 3 

^-C{x\ + x 2 , x\x 2 + xixl) + ^-C(x\x\,x\x 2 + xia;|)) + C(e 2 ) 
9 3 

According to l|12.4fl our reduced 1-dimensional operator has the symbol 

h N ° 

Q E = q e + (D(h+—). 

Put 

q s = sQ e + (1 - s)(q). 
If we assume that e ^> h, then we get according to (|8.61|> : 

£,dx- [ £dx = -ie[ [f{r,a,b,c)} ( p xW - m) dt + 0(e 2 + h), (13.53) 

7i(Q«) ^7l«9» " / 7l«9>) 

where we recall that /^(/^v) £dx is the (real) action along a loop in (q) = Const = (q)(p c ) starting 
and ending at the saddle point p Cl and that J^/q ) £,dx is the corresponding perturbed action for 
Q t . From l|13.46|) . we see that C(xf + minimal precisely on the horizontal circle 

p = 1/2. In the cases B±, E±, D the saddle points belong to {C/,Cb} situated on that circle. 
Since (q) is invariant under reflection in that circle, either the loop 71 ((<?)) is entirely in the upper 
or lower hemisphere intersecting the equator only at p c (and this is happens in the cases B± and for 
one of the saddles in case D) or 71 ((g)) intersects the equator at one more point and is symmetric 
around the equator (and this happens in the cases E± and for one of the saddles in case D). In 
both cases we see that 

/ [C(xj + x 4 ,x 4 1 +x 4 )} ( p x(t) ^dt>0. (13.54) 

Taking into account the form of / in 113.5211 we conclude that for every r the integral in the left 
hand side of l|13.53fl is 7^ except for (a,b,c) in a set of measure 0. For (a,b,c) outside that 
exceptional set, we conclude from the discussion at the end of Section |S] that the spectrum of the 
one dimensional localized operators has a genuinely two-dimensional structure. 



A Proof of Proposition 16.21 

To get a complete normal form we shall do further conjugations with analytic pseudodifferential 
operators of order in such a way that the complete symbol also becomes a function of r, e, h 2 /e 
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and x£. Moreover, we need to do so with errors that are 

0( e -i/(CTO) (rather than merely 
as in ^21- Q is n °t a classical analytic symbol but it has a holomorphic realization and becomes 
a classical analytic symbol, if we allow some of the /i-dependence to appear as an independent 
parameter in the coefficients of the /i-asymptotic expansion. Thus, our starting point will be a 
symbol of the form 

Q = Qo(t, atf, e, h 2 /e) + hQ^r, x, £, e, h 2 /e; h), (A.l) 

where Qi is holomorphic and 0(1) in some fixed complex neighborhood of r — 0, x = £ = 0. 

We define the p-quasi-norm as above, but now it is important that we work in the Weyl 
quantization. To an analytic symbol a we then associate the infinite order differential operator A = 
Op a (x,£,D Xj £]h) as in 13. 5|) . From the definition, we verify the following metaplectic invariance 
property: If K : C 2 ™ — * C 2n is an affine linear canonical transformation and k* , denote the 
usual operation of pull-back and push- forward of functions on C 2 ™, then 

n*Op a K* = Op Kta . (A.2) 

This implies that if we define our quasi-norms with the help of a family of opens sets fit which are 
invariant under k, then 

|||k*o||Ip = \\\o\\\p- (A-3) 

In the case of (|A.lfl . we shall let f2 t be of the form \x\ 2 + |^| 2 < r(t) for a suitable r(t), and we 
observe that these balls are invariant under expisH x ^ when s is real. After applying the inverse 
function of Qq(t, •, e, h 2 /e) to our operator, we may assume that the principal symbol of Q is x£, 
so l|A.l(l simplifies to 

Q = x£ + hQ 1 (T,x,£,e,h 2 /e;h) (A A) 

with a new Q\ having the same properties as the previous one. 

Using the same letters for operators and their symbols, we let Qo = ^(xhD + hDx) be the 
quantization of x£. Notice that exp (2irQo/h) = —I, so that exp (2ir adQ //i) = 1. If B is an 
analytic /i-pseudodifferential operator of order 0, we put 

(B) = _L^ e tQ °/ h Be- tQo/h dt, (A.5) 

and notice that on the symbol level, 

1 [ 2 * 

(B) = — Bo exp UH x fdt. (A.6) 
27r J a 

Also notice that [Qo, (B)] — 0. Choosing the p-quasi-norms as above, we further have that 

\\\(B)\\\ P <\\\B\\\ P . (A.l) 

The equation, 

a,d Qo A = B-(B) (A.8) 

has the solution 

A=\l k{t)e tQ °/ h Be- tQo/h dt = i / k{t)e tQo/h {B - (B)) e - tQ °/ h dt, (A.9) 
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where k(t) is the function with support in [— 7r,7r] which is afhne on [— ir,0[, ]0, 7r] with fc(±7r) = 0, 
fc(±0) = t\- We have 

\\ P <C\\\\{B-{B))\\\ p . (A.10) 



As in Section |3 we see that the map 

A^Ad A (Q) (A.ll) 

has the differential 

5A^ad SA Qo + K(A,SA), (A.12) 

where 

\\\K(A,SA)\\\ p <Cp(h + \\\A\\\ p )\\\SA\\\ p , (A.13) 

under the assumption that = 0(1). 

Now return to i|A.4fl . After a first conjugation, we may reduce ourselves to the case when 
Qi - (Qi) is OQi), so that 

IIIQi - <Qi)|||p < Ch, (A.14) 

for some C > 0, when p < p > 0. We look for A such that Ad^Q commutes with Q , and we try 
A = X)o° with convergence in some p-quasi-norm. Start by solving 

[A , Qo] + hQ x = h(Qi), with \\\A \\\ P < C\\\Q 1 - (Q^W, < O(h). (A.15) 

From l|A~T2|) . l|A~T3)l . we get 

M Ao Q = Q + h(Q 1 ) + hQ 2l (A.16) 

IIIQalllp < h- 1 C P {h + C\\\Q 1 -{Q 1 )\\\ p )\\\Q 1 -{Q 1 )\\\ p (A.17) 

< Cp(l + C 2 )|||Q 1 -(Q 1 )|||, 

< \\\\Qi-(Qi)\\\p, 

where the last estimate holds for < p < po, with po > small enough. Choose Ai with 

[A u Qo] + hQ 2 = h{Q 2 ), pilUp < C\\\Q 2 - (Q 2 )\\\ p < ^|||Qi - <Qi)||| P . (A.18) 

Then 

Ad Ao+Al Q = Q + h{Qi) + h(Q 2 ) + hQ 3 , 
with IIIQslllp < 2 _2 |||Qi — (Qi)|||p- Iterating the procedure, we get Aj with 

\\\Aj\\\ p < C2- j \\\Q 1 - (Qi)\\\p, 

such that if A = Aj'j tncn 

Ad A Q = Q + HQ,) + h(Q 2 ) + HI <^> ||U < 2-ilQx - <Qi>||| p . (A.19) 

The previous discussion shows how to find U so that modulo an error 0{e~ 1 ^ ch ^), U~ 1 QU 
commutes with xhD x . Moreover U~ 1 QU is a classical analytic pseudodifferential operator (after 
allowing h as an independent parameter in the coefficients in the asymptotic expansions). Put 



Ml 



x = e s and work near x = r for some fixed small r > 0. Then xhD x = hD s and since the class 
of analytic pseudodifferential operators is conserved under analytic changes of variables, we know 
that 

U^QU = K„[r,s,hD s ;h), 
where K is an analytic symbol. But [K, hD s ] = 0, so 

K = K..(t, hD s ;h) 

and returning to the x-coordinates, we get the representation □ 



B Study of r j; k 



For simplicity, we restrict the attention to the right half-plane, Re \i > and pick one of the 
equations in (|8.12(l . that we write 

(Im M )ln^-=f(/i), (B.l) 
IMI 

where F (fj,) is uniformly Lipschitz in a neighborhood of 0. As we have already observed, 

1 1 ( \vn 

(iTO^ln— ) = ln^- -j-f- >1, (B.2) 

H ImI V H / 

so (IB. Ill determines a curve of the form 



Imn = /(Re M ), where /'(Re /x) = f - 1 ) « 1. (B.3) 

Vlnl/|(Re/i,/(Re^))|y 

We want to express / in terms of F^Re /u) up to small errors. 

Let us first compare the solution /.i of (|B.1|) with the solution Jl of the simplified equation 



Im/iln — = F(Re u), with Re /i = Re /i. (B.4) 

ImI 

Using that F(fi) — F(Re fi) = 0(lm /it) together with (|B.2f) . we see that 

Im fj, — Im /i — 0( 3—), so Ito/i ~ Im /i. In- — - ~ In -rzrr (B-5) 

top ImI N- 



With this estimate in mind, we now concentrate on the simplified equation i|B.4fl . and we drop the 
tildes for simplicity 

Assume first that we are in the region 

\Im[i\ < 0(Re/x). (B.6) 

Then 



In i n M+0 (») 



Ml V 3 -/ V \i lal/i 
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where we write pL — x + iy. Thus, if [i = fj, solves (IB.4|) and HB.6(I holds, then we first see that 

V ln±' 

X 

and then that 



In 

So, if we assume 



r( 1 + ° (1)( iS> s i^)' (■"> 



LF(a:)| < 0(l)xln-. (B.S 

a; 



then we are in the region (|B.6|) , and the solution /i — /i = x + iy of (|B.4fl takes the form (|B.7(I . 
Combining with (|B.5|> we get under the assumption (|B.8|) . 

(B.9) 



/ 0(1)\F(!) 
/(X) H 1 + lnTJ^nT- 



We next consider the region x <C \y\ <C 1 and assume for simplicity that we have y > 0. Then, 

b R = ta( K 1 + ° (1)< ^i4)' (R10) 



and (IB.4I) takes the form 



y(ln + 0(l)(-f-\ ) = F(x). (B.ll) 

y \ y in-; 



Consider first the simplified problem 



y\n- = z. (B.12) 



With Y = ln(l/y), Z = ln(l/z) (both > 1) we get 

Y-\tlY = Z (B.13) 
Try the approximate solution Yq = Z + In Z. Then by a simple calculation, 

y o -iny o = z + 0(^). 



Since the derivative of the left hand side in IB. 13(1 is close to 1, we see that the solution Y of that 
equation is of the form Y = Y Q + C((ln Z)/Z)\ 

Y = Z+ (l + 0(^y\lnZ. (B.14) 



Hence the solution of (IB. 1211 is of the form 



mini. 



» = e+°(l^r> si- (b.i5) 
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If we replace z by F(x), we get the order of magnitude of the solution to (|B.11|I 



and the assumption that i<j/ reads: 

F(x) 



In „} \ 

F(x) 



The earlier arguments show that this condition is equivalent to 



which indeed is complementary to 1|B.8|1 . Using l)B.16|) . we get 

ln^± mi in — 



t 1 a; In -prr 1 £ In , i 

(*) 2 -^ < o ( i)( J^—L^ < i , 

2/ lni F(ar) ln^v 



where we notice that 

xln -I 



F «1, 



by p,17|) . Hence ljB.ll|l gives 



and applying (|B.15|I with z equal to the right hand side of l|B.19|l . we get 

TJf) 



- n ^^ lnlnl/F , (l + g(l)( £i ^) 2 I^)^ ) 



(B.16) 



> x. (B.17) 



F(x)^>xln-, (B.18) 



ylni = (1 + 0(1)( ^^ ) 2 — ^— )Fix), (B.19) 
y 1111/^(1) 



which simplifies to 

Recall that here Jl = x + iy in the simplified equation. To get the corresponding result for (|B.l 



^ N lnlnlAF\ F(z) 



we apply (|B. 1|) and conclude that (|B.20|) holds for /i = x + iy solving IjB.lf) . under the equivalent 
conditions (|B.17|) . IB. 181) . and assuming also < F(x) < 1, < i « 1. 
Summing up, we have proved: 

Proposition B.l Let F be a uniformly Lipschitz function with \F\ <C 1, defined in a neighborhood 
of £ C. Let /j, = x + if(x) be the solution of Then for small x, we have 

\f'{x)\<0(l)/Ml/\x + if(x)\). (B.21) 

Further, 

f{x) = (1+ h^M ) ln^T' IF(X)1 ~ ° {lM ln(1/|2;|) ' (R22) 

fix) = (1 + 0{1) ^jfp^ ) H f/iF\x)\) ' ^ N ln(1/NI) ' (B ' 23) 
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